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Pressure  sensitive  paints  (PSP)  have  been  used  extensively  in  the  past  decade  to 
measure  surface  pressures  in  steady  flow  conditions.  There  are  many  unsteady 
phenomena  to  which  the  aerospace  community  would  like  to  apply  PSP  as  a 
measurement  technique.  However,  to  date  there  has  not  been  a  viable  model  to  describe 
the  dynamic  behavior  of  the  PSP  system.  A  rigorous  development  fi-om  a  systems 
dynamics  perspective  is  needed  to  advance  PSP  fi-om  a  static  measurement  technique  to  a 
dynamic  one. 

With  this  aim,  two  models  for  the  dynamic  behavior  of  pressure  sensitive  paints 
will  be  presented.  The  first  of  the  two  models  is  a  purely  empirical  approach  to  designing 
a  model  and  compensator.  The  second  model  encompasses  the  physics  of  the  process  by 
which  an  unsteady  pressure  field  over  the  paint  layer  affects  the  layer  and  causes  an 
intensity  of  fluorescence  that  is  fluctuating  in  time.  Within  this  second  model,  two 
different  forms  for  the  static  calibration  are  chosen.  The  first  results  in  purely  linear 
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system  dynamics  whereas  the  second  yields  a  nonlinear  system.  These  models  are  then 
compared  to  experimental  data.  Next,  a  dynamic  compensation  scheme  is  developed  to 
correct  for  the  time-lag  and  amplitude-damping  behavior  of  the  paint.  Finally,  the 
compensation  scheme  is  applied  to  PSP  data  taken  surrounding  an  unsteady  (periodic)  jet 
injected  transversely  into  a  Mach  1.56  freestream. 
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CHAPTER  1 
INTRODUCTION 

Pressure  Sensitive  Paint  (PSP)  has  been  an  emergent  technology  over  the  past  five 
to  ten  years.  Many  different  groups  have  made  contributions  as  is  summarized  in  the 
reviews  by  Bell  et  al.,'  Liu  et  al.,^  and  Lu  and  Winnik.'  Much  success  has  been  achieved 
in  the  application  of  PSP  to  steady  pressure  fields.  In  comparison,  there  are  but  a  few 
papers  discussing  the  application  of  PSP  to  relatively  unsteady  pressure  fields.  It  is  the 
goal  of  this  dissertation  to  present  a  physically  realistic  model  for  the  dynamic  response 
of  PSPs  and  provide  a  dynamic  compensation  scheme  for  PSP  data  analysis.  With  this 
done,  a  dynamic  PSP  experiment  will  be  presented  and  analyzed  using  the  previously 
derived  dynamic  model  and  compensation  scheme. 

Basic  PSP  Phvsics  ■  ' 

PSP  is  a  measurement  technique  that  has  its  roots  in  photochemistry.  It  is  based 
on  the  luminescence  quenching  effect  of  oxygen  on  luminescent  molecules  in  a  polymer 
binder.  Shown  in  Figure  1.1  is  a  schematic  of  a  typical  PSP  layer.  The  schematic  shown 
here  is  for  a  PSP  applied  directly  to  the  substrate  without  a  primer  layer.  The  reason  for 
this  is  that  primer  layers  will  slow  down  the  dynamic  response  of  the  PSP  due  to 
luminophore  molecules  diffusing  into  the  less  permeable  primer.  Because  the  primer 
generally  has  a  lower  oxygen  diffiisivity  than  the  polymer  binder,  the  luminophores 
trapped  in  the  primer  will  exhibit  very  slow  dynamical  responses.  The  photophysical 
process  begins  with  the  excitation.  For  most  PSPs,  the  excitation  is  a  light  source  (strobe. 
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Figure  1.1 --Schematic  of  a  typical  PSP  layer. 

halogen,  LED  or  laser)  with  high  intensity  in  the  blue  to  ultraviolet  region  of  the 
spectrum.  A  given  luminophore  in  its  ground  state  can  absorb  a  photon  and  transition  to 
an  excited  state.'*  Once  in  the  excited  state,  the  luminophore  has  primarily  three 
competing  routes  to  return  to  the  ground  state.  ^  The  first  is  referred  to  as  "radiative 
decay"  or  "luminescence"  which  encompasses  two  phenomena  -  fluorescence  and 
phosphorescence.  For  ruthenium-based  paints,  fluorescence  is  the  primary  radiative 
decay  route  whereas  for  the  platinum  derivative  paints  it  is  phosphorescence.  In  both 
fluorescence  and  phosphorescence,  the  luminophore  molecule  emits  a  photon  and  drops 
back  down  to  the  ground  state.  This  photon  is  generally  in  the  orange  to  red  region  of  the 
spectrum.  The  intensity  of  this  emitted  light  is  the  measurement  variable  of  PSP.  The 
second  decay  route  is  called  "non-radiative  decay"  in  which  the  luminophore  releases  its 
energy  to  the  surrounding  polymer  matrix  in  the  form  of  thermal  energy.  The  third  route 
is  the  critical  one  that  differentiates  an  ordinary  fluorescent  paint  from  a  pressure 
sensitive  paint.  In  this  route,  molecular  oxygen  that  has  permeated  the  layer  may  collide 


3 

with  an  excited  luminophore  and  absorb  its  energy.  This  third  route  is  termed  "oxygen 
quenching." 

Each  of  the  three  methods  of  energy  release  has  an  accompanying  rate  constant  - 
kr,  knr,  kqno2  respectively.  The  first  two  rate  constants,  kr  and  k„r,  have  units  of  1/s .  The 
third  rate  constant  is  composed  of  two  parts,  kq,  which  is  a  constant  with  units  of 
mY(mol  •  s),  and  no2  which  is  the  concentration  of  molecular  oxygen  in  mol/m^  within 
the  PSP  layer.  The  quantum  efficiency  of  PSP  luminescence  is  given  generally  by 

The  quantum  efficiency  is  proportional  to  the  observed  PSP  intensity,  (|)  =  c  •  I ,  with  c 
being  a  proportionality  constant.  If  a  ratio  is  made  of  the  quantum  efficiency  in  the 
absence  of  oxygen  (n^^  =  O)  to  that  in  the  presence  of  oxygen  (no2    O),  the  following 
relation  is  obtained: 


(b      cl  I 

To    _  *^^o  _ 

(|)  ~  cl  ~  I 


--^  =  l  +  K^no2  (1.2a) 


where 

Equation  (1.2a,b)  is  familiarly  known  as  the  linear  Stem-Volmer  equation.  The 
equilibrium  oxygen  concentration  within  the  layer  is  related  to  the  air  pressure  above  the 
layer  via  Henry's  law^ 

"02  =f(Po2)  =  cPo2  =cXo2P  =  tJP-  (1.3) 
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In  the  above  equation,  the  Hnear  form  for  Henry's  law  has  been  presented.  The  solubihty 
of  gases  in  polymers  is  often  modeled  in  this  fashion.''  Substituting  equation  (1.3)  into 
equation  (1.2a)  yields  a  relation  between  the  luminescence  intensity  and  the  air  pressure 
over  the  PSP  layer, 

^  =  1  +  ctK,P.  (1.4) 

In  aerodynamic  testing  it  is  often  not  possible  to  measure  the  luminescent  intensity  in  the 
absence  of  oxygen,  lo.  However,  the  intensity  at  some  reference  pressure,  usually  taken 
to  be  atmospheric,  can  be  measured.  Applying  equation  (1.4)  at  a  reference  condition, 
Iref,  and  at  a  general  condition,  I,  and  ratioing  gives  the  "aerodynamic  testing"  form  of  the 
linear  Stem-Volmer  equation: 

I  ,      1  +  ctK  P  p_o 

^  =  ^  =  J_^.  (1.5) 

I      l  +  aK,P„,  P„,-Q 

In  this  equation,  Q  =  -  l/crK^  is  the  virtual  pressure  (negative  absolute  pressure)  at 
which  the  luminescent  intensity  would  become  infinite,  i.e.,  the  singular  point  of  the 
intensity  fiinction.  The  parameter  Q  is  a  purely  mathematical  construct,  but  has  an 
important  purpose  in  relation  to  the  use  of  equation  (1.5).  Often  equation  (1.5)  is  quoted 
in  the  form''^ 

^  =  a  +  b^.  (1.6) 

However,  in  this  form,  both  the  slope,  b,  and  the  intercept,  a,  will  be  fianctions  of  Pref.  In 
the  form  of  equation  (1.5),  the  parameter  Q  is  independent  of  the  reference  pressure  and 
is  hence  a  more  usable  form.  Typical  values  for  Q  are  in  the  range  of  -  0.1  atm  to 
-  0.2  atm .  A  sensitive  PSP  with  a  steeply  sloped  Stem-Volmer  response  will  have  a  Q 


value  closer  to  zero,  whereas  a  less  sensitive  PSP  will  have  a  more  negative  value  of  Q. 
Solving  for  P  in  equation  (1.5)  yields  the  calibration  for  a  given  PSP  as 

P  =  Q  +  (Pref-Q)Y- 

So,  by  taking  two  measurements,  one  at  a  known  reference  pressure  and  another  at  an 
unknown  pressure  of  interest,  the  intensities  of  the  two  measurements  can  be  transformed 
into  a  single  measurement  of  pressure. 

Today's  modem  CCD  cameras  have  on  the  order  of  a  million  pixels.  One  can 
then  achieve  on  the  order  of  a  million  simultaneous  pressure  measurements  distributed 
over  a  given  model  instrumented  with  PSP  in  a  single  image.  The  physical  spacing 
between  these  measurements  is  essentially  limited  only  by  the  magnification  of  the 
camera  lens.  Put  another  way,  to  achieve  a  denser  set  of  measurements,  one  has  only  to 
zoom  in  tighter  on  a  given  region  of  interest.  Using  microscopes,  sub-micron  resolution 
can  even  be  achieved. 

PSP  in  Steadv  Pressure  Fields      ;    •  .  ; 

Many  experiments  have  been  performed  using  PSP  to  measure  pressures  imder 
steady  conditions.  As  such,  there  are  several  papers  that  review  the  existing  literature.''^'' 
To  summarize  these  reports,  PSP  has  been  used  to  measure  standard  aerodynamic  loads 
on  aircraft  in  both  subsonic  and  supersonic  regimes,  turbomachinery  such  as  fans  and 
compressors,  automobile  wind  tunnel  testing,  as  well  as  numerous  scientific  experiments 
to  investigate  particular  flowfields  of  interest.  One  representative  study  is  that  of  Everett 
et  al.^ 

In  the  paper  of  Everett  et  al.,  a  sonic  jet  of  air  was  injected  transversely  into  a 
supersonic  freestream  at  a  Mach  number  of  1.6.  The  PSP  was  applied  to  the  wall 


surrounding  the  jet  orifice.  The  full-field  wall  pressures  around  the  jet  were  measured  at 
three  different  injectant  pressures.  The  jet-to-crossflow  momentum  flux  ratio, 


took  on  the  values  of  J  =  1.2,  1.7,  and  2.2  for  the  three  injectant  pressures.  In  addition  to 
the  PSP  measurements,  shadowgraph  and  surface  oil-flow  visualization  techiques  were 
used  to  investigate  the  flow.  The  shadowgraph  images  were  useful  in  identifying  the 
shock  structure  generated  as  the  jet  altered  the  flow.  As  the  momentum  flux  of  the  jet 
was  increased  relative  to  the  fixed  freestream,  the  size  and  strength  of  the  shock  structure 
were  seen  to  increase  as  one  would  intuitively  imagine.  With  the  oil-flow,  the 
streamlines  along  the  wall  were  characterized.  From  these  streamlines,  the  locations  of 
the  various  shocks  near  the  wall  could  be  determined.  Placing  static  pressure  taps  in 
strategic  locations  allowed  simultaneous  PSP  and  conventional  pressure  transducer 
measurements.  Using  these  simultaneous  measurements,  an  in-situ  calibration  of  the  PSP 
was  performed.  The  calibration  fiinction  chosen  related  the  indicated  PSP  pressure  to  a 
quadratic  function  in  terms  of  Iref/I.  Spectacular  images  of  the  surface  pressure  were 
achieved  with  a  spatial  resolution  of  85.  l|am/pixel.  Integrating  the  pressure  around  the 
periphery  of  the  jet  gave  a  measure  of  the  effective  back  pressure 


With  increasing  J,  the  effective  back  pressure  was  also  found  to  rise.  Because  of  this,  it 
was  concluded  that  the  effective  back  pressure  was  not  a  constant  ratio  of  either  the  static 
or  stagnation  pressure  behind  a  normal  shock  in  the  freestream  as  was  reported  by  other 
researchers. 


J  = 


(1.8) 


(1.9) 


PSP  in  Unsteady  Pressure  Fields 
The  earliest  reference  of  PSP  in  unsteady  situations  is  a  paper  by  Cox  and  Dunn 
and  contains  the  beginnings  of  a  good  model  for  the  response  of  PSP  to  unsteady  pressure 
fields.'  In  their  paper  they  investigated  oxygen  transport  within  poly(dimethyl  siloxane) 
(PDMS)  via  its  effect  on  the  fluorescence  of  the  luminophore  9,10-diphenyl  anthracene 
(9,10-D).  The  purpose  of  their  work  was  to  determine  such  physical  quantities  as  the 
solubility,  diffusivity  and  activation  energies  of  the  polymer  PDMS.  This  was 
accomplished  with  a  cuvette  of  PDMS  (thickness,  a  =  1cm)  into  which  was  mixed  the 
luminophore  9,10-D.  The  schematic  shown  in  Figure  1.1  is  similar  to  the  geometry  of  the 
cuvette  used  in  their  experiments.  The  cuvette  was  allowed  to  come  to  equilibrium  with 
pure  nitrogen  at  1  atm.  Then,  the  nitrogen  environment  was  replaced  with  pure  oxygen  at 
latm.  Over  the  next  24  hours  the  fluorescence  intensity  at  the  center  of  the  cuvette  (x  = 
4.5mm  to  5.5mm)  was  monitored.  For  PSPs,  the  thickness  is  on  the  order  of  I0[xm  as 
compared  to  the  1cm  thickness  of  the  cuvette.  The  response  time  of  a  PSP  will  be  shown 
to  be  proportional  to  the  thickness  squared.  Thus,  PSPs  work  on  much  shorter  time 
scales  than  in  ref  9.  The  fluorescence  intensity  within  the  cuvette  was  seen  to  fall  off 
quasi-exponentially.  With  the  use  of  a  static  calibration  they  were  able  to  extract  the 
indicated  oxygen  concentration  as  a  fimction  of  time  within  the  PDMS.  They  then 
developed  an  analytical  model  for  the  indicated  oxygen  concentration  using  the  one 
dimensional  diffusion  equation.  For  a  step  change  in  oxygen  concentration  at  the  upper 
surface  (x  =  a)  the  analytical  solution  for  the  oxygen  concentration  within  the  cuvette  is 


n(x,t)-no^^    4y^(  1 

'        "  i=0  ^ 
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This  model  was  found  to  accurately  represent  the  experimental  data.  The  regression 

analysis  yielded  values  of  diffusivity  between  D  =  2.0x  10"'cm^  /s  at  T  =  5°C  and 

D  =  6.1x  lO'^cm^  /s  at  T  =  45°C  for  oxygen  in  PDMS.  However,  in  the  general  use  of 
PSPs  the  coating  is  viewed  from  above  rather  than  from  the  side  as  was  done  in  this 
paper.  That  is,  the  intensity  measured  from  a  PS?  is  the  integrated  intensity  of  the  entire 
layer.  The  integrated  intensity  as  a  fiinction  of  time  for  an  arbitrary  unsteady  pressure 
field  will  be  shown  later  to  have  a  troublesome  nonlinearity. 

The  work  of  Mills  and  Chang'"  was  the  first  to  look  at  the  dynamic  response  of 
optical  sensor  films  (to  be  later  called  PSPs).  The  response  and  recovery  times  of  the 
optical  signal  (luminescence  output)  from  the  coating  were  measured.  The  response  time 
was  defined  as  the  time  required  for  the  optical  signal  to  equilibrate  to  a  positive  step  in 
oxygen  concentration,  and  the  recovery  time  was  the  time  required  for  equilibrium  due  to 
a  negative  step  in  concentration.  It  was  shown  that  the  difference  in  the  response  and 
recovery  times  was  related  to  the  hyperbolic  fiinction  between  the  fluorescence  intensity 
and  oxygen  concentration.  This  paper  was  the  first  attempt  at  characterizing  the 
dynamics  of  PSP  and  used  a  time-domain-based  technique. 

A  paper  similar  to  that  of  Mills  and  Chang'"  was  that  of  Baron  et  al."  Here,  the 
inverse  of  the  PSP  intensity  was  compared  to  the  pressure  as  the  pressure  was  quickly 
stepped  from  near  vacuum  up  to  latm  of  standard  air.  In  this  work  the  assumption  was 
made  that  the  PSP  responded  as  a  first  order  dynamical  system  with  up  to  two  time 
constants.  The  given  dynamical  system  was  numerically  convolved  with  the  pressure 
signal  from  a  fast  conventional  fransducer  to  yield  a  pressure  that  was  compared  to  the 
pressure  indicated  by  the  PSP.  Using  nonlinear  least-squares  curve  fitting  the  time 


constants  and  term  weights  were  determined  yielding  an  empirical  model  for  the  dynamic 
response  of  the  PSPs.  With  these  empirical  models,  they  were  able  to  get  extremely  good 
comparisons  between  curvefits  and  experimental  data. 

The  work  of  Borovoy  et  al.'^  was  the  first  actual  wind  tunnel  application  of  PSP 
in  quasi-unsteady  conditions  found  in  the  literature.  In  this  experiment  the  pressure  on  a 
cylinder  in  Mach  6  flow  was  investigated.  To  achieve  the  Mach  6  flow,  the  cylinder  was 
located  perpendicular  to  the  exit  of  a  shock  tube.  It  is  stated  that  the  99%  relaxation  time 
of  a  PSP's  fluorescence  is  given  by  the  equation 

12  a^ 

to.99=— —  •  (1.11) 
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Once  the  shock  had  passed  the  cylinder,  the  timing  system  in  the  experiment  would  wait 
an  amount  of  time  given  in  equation  (1 . 11)  and  then  trigger  a  single  flash  from  a  xenon 
flash  lamp.  The  PSP  data  were  compared  to  an  analytical  solution  of  Euler's  equations 
for  nonviscous  flow  and  found  to  agree  to  within  10%.  It  should  be  noted  that  this 
application  was  not  truly  an  example  of  unsteady  PSP,  as  the  data  were  taken  after  the 
shock  had  already  passed  the  cylinder  and  the  PSP  had  been  given  sufficient  time  to 
respond  to  the  step  change  in  pressure.  That  is,  the  flow  field  was  unsteady  in  time  but 
not  during  the  period  over  which  the  data  were  acquired. 

In  separate,  but  concurrent  developments,  the  response  of  PSP  to  periodic 
pressure  fields  was  investigated  by  Engler'^  and  Carroll  et  al.'''  In  Engler's  work,  the  PSP 
was  subjected  to  periodic  forcing  at  frequencies  from  0.1  Hz  to  50Hz.  The  main  goal  of 
this  investigation  was  the  characterization  of  the  dynamic  range  and  the  pressure 
resolution  of  the  PSP.  No  effort  was  made  here  to  perform  dynamic  compensation  or  to 
try  to  mathematically  or  physically  define  the  dynamics  of  the  PSP  system.  Carroll  et  al."* 


presented  experiments  similar  in  form  to  those  of  Engler,  namely,  the  response  of  PSP  to 
sinusoidal  pressure  fields.  Much  of  the  analysis  performed  here  was  done  in  the 
frequency-domain.  It  was  shown  that  the  PSP  responds  as  approximately  a  first  order 
dynamic  system  based  on  the  amplitude  response.  The  phase  response  of  the  PSP  system 
did  not  match  well  to  the  standard  first  order  system  model. 

Carroll  et  al.'^  investigated  the  response  of  PSP  to  a  step  change  in  pressure.  In 
this  paper  are  the  beginnings  of  a  physically  realistic  model  for  the  dynamic  response  of 
PSP.  In  summary,  the  paper  reports  on  an  analytical  model  based  on  the  one-dimensional 
mass  diffusion  equation  to  describe  the  fransport  of  oxygen  within  the  PSP  layer.  This 
model  was  nonlinear  and  was  implemented  numerically.  In  addition,  empirical  models 
based  on  a  single-term  first  order  response  and  on  a  two-term  first  order  response  were 
developed.  Using  nonlinear  least-squares  curvefitting  the  two  empirical  models  were  fit 
to  the  data.  The  single-term  first  order  model  was  shown  to  inadequately  fit  the  data. 
The  two-term  and  diffusion-based  models  were  both  shown  to  have  excellent  agreement 
with  the  experimental  data.  An  equivalent  form  of  the  diffusion-based  model  will  be 
derived  in  its  entirety  in  a  later  section.  *  ! 

Winslow  et  al.'^  developed  a  linearized  dynamics  model  for  PSP  based  on  the 
assumption  that  the  luminescence  intensity  from  the  paint  is  proportional  to  the  average 
oxygen  concentration  within  the  layer.  The  frequency  response  of  this  model  is  that  of  a 
"1/2-order"  dynamic  system  with  an  amplitude  response  of-lOdB/decade  above  the 
cutoff  frequency  and  a  high  frequency  phase  shift  of -45°.  Curvefitting  this  model  to 
PSP  data  in  the  frequency  domain  gave  a  value  for  the  time  constant  of  the  layer.  A 
dynamic  compensator  was  then  developed  by  inverting  the  frequency  domain  response 
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and  performing  an  inverse  Fourier  transform.  The  compensator  took  on  the  form  of  a  six 
term  finite  impulse  response  (FIR)  filter.  Applying  this  compensator  to  experimental 
data  yielded  an  impressive  increase  in  the  dynamic  response  of  the  PSP.  This  was  the 
first  instance  of  dynamic  compensation  of  PSP  data. 

Nearly  concurrent  to  the  paper  of  Carroll  et  al.'^  was  that  of  Masoumi  et  al." 
Effectively,  the  same  dynamical  model  was  used  in  both  works.  In  the  latter  work,  the 
luminescence  response  of  PSPs  to  positive  and  negative  steps  in  pressure  was  examined. 
The  coating  thicknesses  were  on  the  order  of  300|im.  Upon  performing  a  curvefit  of  their 
dynamics  model  to  the  data,  values  for  the  diffusivity  of  1.0  x  10"^  cm  Vs  for  platinum 
octaethylporphine  (Pt(OEP))  in  poly(n-butylaminothionylphosphazene)  (PATP)  and 
1.4x10"'  cm  /  s  for  ruthenium  biphenyl  phenanthroline  chloride  (Ru(ph2phen)3Cl2)  in 
PATP  were  obtained. 

Carroll  et  al.'^  reapplied  the  dynamics  model  of  Winslow  et  al.'^  to  investigate  the 
diffusivity  of  PSP  coatings  as  a  function  of  temperature.  A  PSP  coating  was  subjected  to 
step  increases  in  pressure  over  a  range  of  temperatures.  From  the  time  constants 
measured  at  the  different  temperatures,  values  for  the  mass  diffusivity  were  calculated. 
Then,  modeling  the  diffusivity  using  an  Arrhenius  relation 


D  =  Dp  exp 


(1.12) 


rt; 

a  value  of  2.48kcal/mole  was  found  for  the  activation  energy,  E,  of  the  coating.  This 
value  compared  well  with  values  fi-om  other  sources. 

Hubner  et  al.'^  performed  PSP  experiments  on  a  flat  plate  and  a  wedge  in  a  shock 
tube.  A  first  series  of  tests  used  a  flat  plate  painted  with  PSP,  a  continuous  illumination 


source,  and  a  PMT  to  measure  the  integrated  light  output  from  the  whole  of  the  plate.  As 
the  initial  shock  passed  the  plate,  the  PSP  was  subjected  to  a  step  in  the  pressure.  When 
the  reflected  shock  came  back  past  the  plate,  the  pressure  stepped  up  even  more.  From 
this  data,  the  model  of  Winslow  et  al.'^  was  used  to  calculate  a  time  constant  of  3.52ms 
for  the  coating.  A  second  series  of  tests  involved  a  wedge  obstructing  the  flow.  A 
photographic  flash  lamp  was  used  to  excite  the  paint  and  a  CCD  camera  acquired  a  single 
image  for  each  run  of  the  shock  tube.  The  results  showed  a  pressure  gradient  on  the 
wedge  due  to  the  effective  area  reduction  within  the  duct.  Overall,  it  was  demonsfrated 
that  thin  films  of  PSP  could  be  used  to  measure  pressures  in  short-duration  facilities. 

The  final  paper  to  be  discussed  here  was  from  Hubner  et  al.^°  Pressure  and 
temperature  measurements  were  made  on  an  elliptic  cone  model  in  Mach  7.5  flow.  The 
tests  were  performed  at  the  48-inch  hypersonic  shock  tunnel  at  Calspan-University  of 
Buffalo  Research  Center.  Due  to  the  short  duration  of  the  flow,  a  method  similar  to  that 
of  Borovoy  et  al.'^  was  used  wherein  the  PSP  was  allowed  to  come  to  equilibrium  with 
the  flow  and  then  a  photographic  strobe  was  used  to  illuminate  the  paint.  The  data 
showed  that  the  pressure  on  the  surface  of  the  model  ranged  from  3  to  12  kPa.  A  second 
set  of  tests  was  performed  with  a  high  frame  rate,  8-bit  camera.  In  these  tests  the  camera 
acquired  several  images  of  a  temperature  sensitive  paint  (TSP)  coated  model  during  a 
single  flash  from  the  strobe.  The  temperature  data  from  these  tests  were  used  to  compute 
heat  transfer  rates  on  the  surface  of  the  model.  The  results  compared  well  to 
measurements  made  with  thin-film  heat  transfer  transducers. 

As  can  be  seen  in  the  literature  to-date  there  have  been  many  papers  investigating 
sinusoidal  or  step-like  pressure  fields,  but  none  that  have  presented  a  PSP  dynamics 


model  in  a  fashion  for  use  with  arbitrary  pressure  fields.  It  is  the  goal  then  of  this 
dissertation  to  present  a  sound  mathematical  and  physical  model  for  PSP  dynamics. 
Moreover,  a  compensation  scheme  will  be  designed  based  on  the  dynamics  models  to 
correct  for  the  amplitude-damping  and  phase-shifting  behavior  of  the  paints.  Finally,  a 
ftiU-field  dynamic  experiment  of  PSP  will  be  presented  along  with  the  results  of  the 
compensation  scheme  having  been  applied  to  the  data. 

Transverse  Jets 

The  study  of  a  jet  of  air  injected  transversely  to  a  supersonic  freestream  is  a 
classical  problem  that  has  been  the  subject  of  many  studies.  A  schematic  of  the  general 
flowfield  is  shown  in  Figure  1.2.  The  main  purposes  for  which  this  flowfield  have  been 
studied  is  in  relation  to  thrust  vector  control  and  to  ftiel  injection  ports  in  scramjet 
engines.  For  thrust  vector  control,  the  integrated  wall  pressure  is  of  primary  importance. 
For  scramjet  engines  it  is  important  for  the  fiiel  (injectant)  to  become  fiilly  mixed  with 
the  freestream  air  in  as  short  a  time/distance  as  possible.  Additionally,  because  the 
injectant  is  to  be  burned  downstream,  it  is  necessary  to  push  the  injectant  as  far  from  the 
wall  area  as  possible  to  reduce  heating  of  the  wall.  Hence,  the  goal  of  many  studies  has 
been  on  determining  the  wall  pressure  distribution  or  on  finding  optimum  parameters  to 
increase  mixing  or  to  maximize  penetration  of  the  jet.  What  follows  is  a  brief 
introduction  to  the  existing  literature  on  the  evolution  of  transverse  jets  in  supersonic 
crossflow. 

The  first  work  to  be  published  on  the  transverse  jet  in  supersonic  crossflow  came 
from  Cubbison  et  al.  out  of  the  NASA  Lewis  Research  Center  in  1960.^^  In  their  paper, 
extensive  wall  pressure  measurements  were  made  surrounding  a  0.062"  jet.  In  this  study 
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Figure  1 .2~Transverse  injection  flowfield  schematic  (Gruber  et  al.^'). 

there  were  basically  three  independent  parameters  which  could  be  varied  -  jet  pressure, 
freestream  Mach  number,  and  pressure  altitude.  The  range  of  jet  pressures  was  from  50 
to  440  psig.  Three  freestream  Mach  numbers  of  2.92, 4.84  and  6.4  were  investigated. 
The  pressure  altitude  was  fiirther  varied  from  55,000  to  1 15,000  feet.  This  range  of 
pressure  altitudes  combined  with  the  range  of  Mach  numbers  resulted  in  a  range  of 
Reynolds  numbers  between  0.84x10^  and  7.78x10^  per  foot.  The  wall  surrounding  the  jet 
was  instrumented  with  on  the  order  of  100  pressure  taps  arrayed  upstream,  downstream 
and  spanwise  relative  to  the  jet  orifice.  The  general  flowfield  resulted  in  a  high  pressure 
region  immediately  in  front  of  the  jet  orifice  and  a  low  pressure  region  behind  the  jet. 
The  results  of  this  study  indicated  that  the  jet-to-freesfream  pressure  ratio  had  a  large 
effect  on  the  pressure  distribution  on  the  surface  surrounding  the  jet.  The  pressure 
distribution  was  also  affected  by  changing  the  freestream  Mach  number  at  a  constant 


15 

pressure  ratio.  The  effect  of  the  Reynolds  number  was  small  compared  to  the  preceding 
two  effects. 

In  a  later  examination  by  Schetz  et  al.   the  location  of  the  center  of  the  first  Mach 
disk  was  correlated  to  a  pressure  ratio  between  the  jet  pressure  and  an  "effective  back 
pressure."  This  effective  back  pressure  was  an  analogy  to  the  back  pressure  seen  by  a  jet 
issuing  into  a  quiescent  medium.  The  value  chosen  for  this  effective  back  pressure  was 
defined  as  80%  of  the  static  pressure  behind  a  normal  shock  in  the  fi-eestream. 

A  paper  by  Orth  et  al.^''  looked  at  the  effect  of  changing  the  shape  of  the  injector 
port  as  well  as  the  angle  of  injection  into  the  freestream.  A  control  volume  analysis  of 
the  fluid  in  the  jet  was  also  presented.  In  this  study,  an  effective  back  pressure  equal  to 
2/3  of  the  fi-eestream  stagnation  pressure  was  used.  It  was  argued  that  the  effective  back 
pressure  should  be  the  average  pressure  field  around  a  cylindrical  obstruction  on  a  flat 
plate.  So,  the  definition  used  here  differed  from  that  of  Schetz  et  al.^''  Several 
experiments  were  performed  to  better  understand  the  factors  affecting  penetration  of  the 
jet  into  the  freestream.  A  first  series  of  tests  used  a  circular-ended  rectangular  slot 
oriented  along  the  flow  direction  and  perpendicular  to  the  flow  as  well  as  an  equal  area 
circular  port  as  the  injector  shape.  These  tests  were  run  at  a  Mach  number  of  2. 1  and  a  jet 
to  back  pressure  ratio  of  5.23.  It  was  found  fi-om  these  tests  that  the  shape  of  the  injector 
port  had  no  real  effect  on  the  penetration  of  the  jet  as  was  originally  hoped.  A  second 
series  of  tests  was  a  parametric  study  of  the  injectant  pressure  and  Mach  number.  It  was 
found  that  the  distance  from  the  wall  of  the  center  of  the  first  Mach  disk  was  related  to 
the  jet-to-freestream  pressure  ratio. 


Billig  et  al.    and  Heister  and  Karagozian  performed  detailed  vortex  modeling 
of  the  jet/crossflow  interaction.  These  models  were  used  to  predict  the  effects  of  varying 
the  freestream  Mach  number,  the  velocity  of  the  jet  as  it  exits  the  wall,  and  the  Mach 
number  of  the  jet  in  relation  to  the  amount  of  penetration.  It  was  determined  that  the 
primary  variable  of  interest  for  enhancing  penetration  was  increasing  the  jet-to-crossflow 
momentum  ratio  (see  equation  1.8). 

As  new  measurement  techniques  were  developed,  new  visual  studies  of  the  jet  in 
crossflow  problem  were  examined.  Among  the  measurement  techniques  employed  are 
planar  laser-induced  fluorescence  (PLIF),  '   classical  Schlieren  images,   planar  Mie 
scattering,^'  and  laser  Doppler  velocimetry  (LDV).^°  Each  of  these  methods  was  used  to 
primarily  to  investigate  penetration  and  spreading  of  the  jet.  The  latest  measurement 
technique  used  to  attack  the  jet  in  crossflow  problem  has  been  PSP  as  was  described  in 
the  earlier  section  "Static  PSP." 

It  has  been  shown  by  the  works  of  Johari  and  Paduano,^'  Hermanson  et  al.,^^  and 
Johari  et  al.   that  both  penetration  and  mixing  can  be  enhanced  through  the  use  of  a 
pulsed  or  modulated  jet  in  crossflow.  It  should  be  noted  that  these  studies  were  all 
performed  in  an  incompressible  medium  as  opposed  to  the  compressible  flows  we  have 
previously  discussed.  However,  the  analogy  may  still  stand.  If  this  is  the  case,  it  should 
be  of  interest  to  investigate  a  pulsed  jet  in  supersonic  crossflow.  This  then  will  be  the 
form  of  the  dynamic  flowfield  to  which  unsteady  PSP  measurements  will  be  tested.  The 
results  of  these  tests  will  be  presented  in  a  later  chapter  of  this  dissertation. 


CHAPTER  2 
PSP  MODEL  DEVELOPMENT 

The  most  important  step  in  performing  a  dynamic  analysis  of  PSP  data  is  model 

selection.  What  follows  are  two  different  approaches  yielding  models  increasing  in 

design  complexity  and  accuracy.  The  first  approach  is  purely  empirical  in  nature.  The 

second  is  based  on  solving  the  one-dimensional  mass  diffusion  equation  in  conjunction 

with  the  relation  between  intensity  and  pressure  under  static  conditions. 

Empirical  Models 

Perhaps  the  easiest  model  that  can  be  developed  is  one  which  does  not  even 
attempt  to  understand  the  physics  of  the  system  in  question,  but  rather,  simply  measures 
the  input  and  output  of  the  system  and  uses  the  data  as  the  characterization.  This 
approach  is  often  referred  to  as  the  "black  box"  approach.  All  one  must  do  is  define  the 
input/output  mapping  for  the  system  as  a  function  of  frequency  in  order  to  create  a  crude 
compensator  to  correct  for  the  amplitude  damping  and  phase  shifting  behavior.  The  only 
underlying  assumptions  to  this  type  of  model  are  that  the  system  being  modeled  must  be 
linear  and  time-invariant.  That  is,  for  an  input  given  by 

Pin(t)  =  Po+P,sin(cot).  (2.1) 
the  output  from  the  system  must  be  of  the  form 

Pou.  W  =  Po  +  P,  A(co)sin(tot  +  (t)(co)).  (2.2) 
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In  the  above  equation  A  and  ^  are  the  ampHtude  ratio  and  phase  shift,  respectively,  of  the 
output  relative  to  the  input  and  are  functions  of  the  driving  frequency  co.  Equations  (2.1) 
and  (2.2)  can  be  more  easily  related  to  one  another  in  the  frequency  domain  via 

Pou.  M  =  HO«)Pi„  M  =  A(co)e-J*('»)p,„  (jco) .  (2.3) 
Shown  in  Figure  2.1  is  the  frequency  response  of  a  16|im  thick  PSP  sample  based  on 
input/output  data.  This  data  was  acquired  using  the  experimental  setup  used  by  Winslow 
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et  al.   and  Carroll  et  al.   The  general  uncertainties  of  the  measurements  are  ±0.1  kPa  for 
the  conventional  pressure  transducer  and  ±2%  of  the  reading  for  the  PSP  measurements. 
Based  on  the  experimental  data,  a  transfer  fiinction  of  the  form 

(2.4) 

1  +  J(0X2 

was  chosen  to  model  the  system.  This  particular  form  of  transfer  fiinction  is  that  of  a  lag 
system.  The  transfer  fiinction  was  applied  to  the  input  data,  Pin(t),  from  a  high  frequency 
response  pressure  fransducer  using  the  Matlab  computer  software  to  achieve  a  model 
output  pressure,  Pmodei(t).  The  parameters  xi  and  T2  were  varied  until  the  root  mean 
square  error  (RMSE)  between  the  PSP-indicated  pressure,  Pout(t),  and  the  model  output 
was  minimized  in  the  time  domain.  The  resulting  values  of  ti  and  X2  were  0.012  and  0.15 
seconds  respectively.  The  model  fransfer  fiinction  is  also  shown  in  Figure  2.1.  Figure 
2.2  contains  the  time  domain  data  for  this  experiment.  Shown  in  the  figure  are  the  PSP- 
indicated  pressure,  the  pressure  measured  by  the  conventional  transducer,  the  model 
output  data  and  the  residual  between  the  PSP  indicated  pressure  and  the  model  output. 
The  absolute  error  between  the  PSP  data  and  the  model  output  data  is  less  than  2  kPa 
(RMSE  =  0.97  kPa).  As  can  be  seen,  the  transfer  fimction  of  equation  (2.4)  does  a 
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reasonably  good  job  of  approximating  the  PSP  system  dynamics.  However,  the 
parameters  of  this  curvefit  are  dependent  upon  the  thickness  and  chemistry  of  the  PSP 
layer  in  an  unknown  manner.  That  is,  if  the  thickness  or  formulation  of  the  PSP  were 
changed,  a  different  transfer  function  would  result.  Hence,  the  empirical  approach  is 
useful  for  a  single  sample,  but  does  not  transfer  well  to  the  type  of  experiments  one  might 
actually  want  to  perform.  Therefore,  a  more  thorough  imderstanding  of  the  physical 
processes  involved  in  the  dynamics  of  PSP  is  required. 

Diffusion  Based  Models 

As  was  mentioned  earlier,  PSP  is  based  on  the  luminescence  quenching  effect  of 
oxygen  on  luminescent  molecules  in  a  polymer  binder.  This  model  tracks  the  process  by 
which  an  unsteady  pressure  above  the  coating  affects  the  oxygen  distribution  within  the 
PSP  and  in  turn  causes  the  intensity  of  luminescence  to  change  in  time. 

Under  static  conditions,  the  oxygen  concentration  within  the  layer  will  be  constant 
and  related  to  the  pressure  over  the  surface  via  a  linear  sorption  law  (Henry's  Law) 

n  =  ctP  .  (2.5) 
Also  under  static  conditions,  the  intensity  of  luminescence  is  related  to  the  pressure  over 
the  surface  via  a  static  calibration 

I  =  f  (P)  (2.6a) 

or 

P  =  g(l)  (2.6b) 
due  to  the  luminescence  quenching  effect  of  oxygen  within  the  layer.  Note  here  that  the 
functions  f  and  g  are  inverse  functions.  For  the  moment  a  generic  relation  between 
pressure  and  intensity  is  assumed.  Later,  specific  forms  for  the  static  calibration  function 
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will  be  used  to  complete  the  system  dynamics.  Letting  l'  =I/Iref  and  P*  =  P/P^f  the 
calibration  equation  can  be  brought  into  a  dimensionless  form 

r=f(p*).  (2.7) 

In  the  above  equation,  the  non-dimensional  integrated  intensity  I*  can  be  replaced  by  the 
dimensionless  intensity  per  unit  depth,  J* ,  by  noting  that 

T'      I       I/a  J 

I  =r-  =  rV  =  -i-  =  J  •  (2.8) 

Also,  pressure  can  be  replaced  with  oxygen  concentration  by  noting  that 

P       aP       n  . 

P  =  =  =  =  n  .  (2.9) 

A  static  calibration  equation  relating  the  intensity  per  unit  depth  to  the  oxygen 
concentration  can  then  be  given  as 

J*=f(n').  (2.10) 
Neglecting  any  internal  light  attenuation,  this  equation  also  holds  on  a  local  level,  that  is, 
at  any  point  within  the  layer  (Cox  and  Dunn^).  Therefore,  equation  (2.10)  can  be  used 
within  the  layer  by  considering  both  J '  and  n*  to  be  functions  of  depth,  x.  In  addition, 
this  equation  holds  as  a  function  of  time  assuming  that  the  kinetics  of  the  luminescence 
and  quenching  processes  are  much  faster  than  the  kinetics  of  the  diffusion  process.  The 
kinetics  of  the  luminescence  and  quenching  processes  are  on  the  order  of  microseconds 
for  the  ruthenium-based  paints  (Bell  et  al.,'  Liu  et  al.,^  Schanze  et  al.,'^  Bums  and 
Sullivan^")  and  on  the  order  of  lO's  to  lOO's  of  microseconds  for  the  platinum-based 
paints  (Bell  et  al.').  Then,  since  the  kinetics  of  the  diffusion  process  are  on  the  order  of 
milliseconds  to  seconds  for  typical  coating  thicknesses  (Baron  et  al.",  Winslow  et  al.'^ 
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Carroll  et  al.'^  Hubner  et  al.'^),  the  assumption  is  indeed  valid.  Restating  equation  (2.10) 
with  the  new  dependencies  one  observes  that  the  static  calibration  provides  a  dynamic 
relation  between  intensity  at  a  given  depth  within  the  PSP  layer  and  the  oxygen 
concentration  at  that  depth, 

J*(x,t)=f(n*(x,t)).  (2.11) 
The  next  step  is  to  solve  for  the  oxygen  concentration  as  a  function  of  depth  and  time  for 
a  given  surface  pressure  input,     (t) .  This  will  be  done  using  the  one-dimensional  mass 
diffusion  equation  with  boundary  conditions  and  initial  conditions  of 

a'n    1  an 

n(a,t)  =  aPjt)  (2.12b) 

an(o,t)  ^ 

-^  =  0  (2.12c) 

n(x,0)  =  aP.„(0).  (2.12d) 
The  governing  equation,  (2.12a),  assumes  that  the  PSP  layer  is  infinite  in  the  transverse 
directions  and  that  the  diffusivity  is  constant.  The  first  boundary  condition,  equation 
(2.12b),  makes  the  assumption  that  the  kinetics  of  the  sorption  process  are  much  faster 
than  the  kinetics  of  the  diffusion  process  so  that  the  surface  of  the  PSP  is  in  equilibrium 
with  the  air  above.  The  second  boundary  condition,  equation  (2.12c),  represents  a  non- 
penetration  condition  at  the  substrate.  The  initial  condition,  equation  (2.1 2d),  states  that 
the  oxygen  concentration  across  the  layer  is  initially  constant  (i.e.  in  equilibrium 
initially).  The  set  of  equations  (2. 12a-d)  or  their  underlying  assumptions  have  been  used 
in  many  other  analyses  of  the  solubility  and  difftision  of  gases  in  polymers.''-'"''^"'^'^^ 
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These  equations  can  be  non-dimensionalized  by  letting  x'  =  x/  a,  t*  =Xlx  where 
T  =  a  V  D ,  and  making  use  of  equation  (2.9)  to  arrive  at 

n*(l,f)=P:(f)  (2.13b) 

5n*(0,t')  ^ 

J^:   ^  =  0  (2.13c) 

n*(x*,0)=P*(0).  (2.13d) 
In  PSP  measurements,  the  luminescence  intensity  is  the  quantity  actually  measured.  Thus 
it  is  instructive  to  recast  equations  (2.13a-d)  in  terms  of  intensity.  Inverting  equation 
(2.1 1)  to  get  n*(x*,t')=  g(j*(x*,t*))  and  substituting  into  equations  (2.13a-d)  yields 


^*  (5x7""  5g/ar 


(2.14a) 


ax  j 

r(i,t-)=f(p:(t-))  (2.14b) 

aj*(o,t*)  , 

-^  =  0  (2.14c) 

J-(x',0)=f(p:(0)).  (2.14d) 
Equations  (2.14a-d)  represent  a  nonlinear  partial  differential  equation  with 
nonhomogeneous  boundary  conditions  describing  the  intensity  distribution  as  a  function 
of  depth  and  time.  Notice  that  recasting  the  formulation  in  terms  of  intensity,  J* ,  instead 
of  oxygen  concentration,  n* ,  introduces  a  nonlinear  term  (last  term  on  right-hand  side  of 
2.14a).  The  problem  can  be  restated  via  a  change  of  variables 

K(x-,t-)=r(x-,t-)-f(p:(t-))  (2.15) 
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to  achieve  a  nonlinear  partial  differential  equation  with  homogeneous  boundary 
conditions  and  initial  conditions 


ydx  J 


at* 


^*   (ax*)"  5g/aj* 

K(l,t*)=0 

ax* 
k(x',o)=o. 

The  output  of  a  PSP  measurement  is  the  integrated  intensity  (neglecting  internal 
absorption) 

i:,(t*)=jj*(x*,t*)dx*=|K(x*,t*}lx*+f(p*(t*)). 


(2.16a) 

(2.16b) 
(2.16c) 
(2.16d) 


(2.17) 


The  integrated  intensity  is  then  related  to  the  indicated  pressure  through  the  static 
calibration,  P:.(t*)=g(i:„,(t*)). 

Equations  (2.16a-d)  are  an  extremely  important  set  of  equations  representing  the 
difference  between  the  intensity  output  by  the  uppermost  layer  of  the  PSP  and  the 
intensity  output  from  a  given  depth  within  the  PSP  layer.  The  term  on  the  left-hand  side 
(LHS)  and  the  first  term  on  the  right-hand  side  (RHS)  of  equation  (2. 16a)  are  both  linear 
terms.  The  last  term  on  the  RHS  of  that  same  equation  is  a  forcing  term.  The  middle 
term  on  the  RHS  is  in  general  highly  nonlinear  and  contains  a  forcing  part  as  well. 
Because  of  the  nonlinearity  of  this  term  it  is  not  possible  in  general  to  apply  classical 
controls  methods  as  will  be  shown  later.  Using  two  different  assumed  calibration 
functions  and  modal  analysis,  a  system  model  now  will  be  developed  and  simulated. 
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Diffusion  Based  Model  with  Linear  Calibration 

The  first  calibration  function  to  be  investigated  will  be 

l'=f(p')=Co+c,P'  (2.18a) 


or 


P*=g(I*)  =  ^-^.  (2.18b) 
c, 

This  form  of  calibration  is  not  typically  used  in  PS?  analysis;  however,  it  has  one 
important  feature.  Using  a  calibration  which  is  linear  between  intensity  and  pressure 
causes  the  (normally)  nonlinear  term  from  equation  (2.16a)  to  drop  out  yielding 

aK(x-,t-)  a'K(x-.t-)  af(p:(t-)) 
St-    ~  (ax7  a- 

which  is  now  a  linear  partial  differential  equation  with  forcing.  The  solution  to  this 
equation  will  be  assumed  via  separation  of  variables  as 

K(x-,t-)=|;a,(t>,(x-)  (2.20) 

i=l 

where     (x* )  are  sine  or  cosine  functions.  This  form  of  solution  is  based  on  a  method 
called  modal  analysis  and  is  often  used  in  linear  systems  dynamics  modeling.  The  spatial 
eigenfiinctions  ^i{x')  must  satisfy  the  boundary  conditions  (2.16b,c)  giving 

4'i(x')=cos(XiX*)  where     =(2i-l)^.  (2.21) 

Substituting  into  equation  (2.19),  multiplying  by  'Pk(x')  and  integrating  across  the 
thickness  of  the  layer  gives 


d,(.-)=-X;a,(t-)-2^)f(p;(,-))  (2.22) 
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Substituting  equations  (2.20)  and  (2.21)  into  equation  (2.17)  yields  the  ou^ut 

C(.-)=Za,(,-)^+f(p-(.-)).  (2.23) 

i=l  '^i 

In  implementation  the  infinite  sum  of  equation  (2.20)  will  be  truncated  after  some 
number  of  terms,  N.  After  truncation,  equations  (2.22)  and  (2.23)  can  be  written  in  the 
standard  linear  systems  form 


/  : 

0 

0' 

f  : 

0 

-^^ 

0 

+ 

V  • 

J 

0 

0 

V  • 

0  -2sm{h,)/X; 


f(p:.(f))' 


(2.24a) 


C(f)= 


sm 


V(p:(.-))J 


(2.24b) 


This  system  is  a  multi-input,  single-output  system  with  inputs  of  f  (p*  (t* ))  and  f  (p*  (t' )), 
state  vector     (t* )  and  output  C  (t* ).  A  useftil  property  of  this  model  is  that  the  system 
matrix  is  diagonal  with  distinct,  nonzero  eigenvalues.  The  usefiilness  of  this  property  is 
m  calculating  the  observability  and  controllability  matrices.  It  can  be  shown  that  these 
matrices  are  both  fiill  rank,  thus  implying  that  the  system  is  both  observable  and 
reachable.  These  two  properties  should  allow  a  simple  compensator  to  be  constructed. 
There  is  one  additional  point  of  interest  concerning  this  system.  Since  it  is  a  linear 
system,  and  the  relation  between  intensity  and  pressure  is  linear,  equations  (2.24a,b)  can 
be  rewritten  via  a  linear  transformation  as  a  system  between  the  input  pressure  and  the 
indicated  output  pressure 


/  : 

0 

0" 

f  : 

•) 

0 

0 

a.(l 

+ 

V  • 

J 

0 

0 

V  • 

) 

0  -2sin(Xi)/Xi 


p-(.-)- 


(2.25a) 


p:.(t*)= 


sin(A,j ) 


f    :  \ 


a,(f) 
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(2.25b) 


It  can  be  shown  that  the  above  output  is  identical  to  the  solution  for  n'(x*,t')  integrated 
across  the  layer,  that  is, 

I 

Po*u.(t*)=  Jn*(x',t')dx*  (2.26) 

0 

where  n*(x',t*)  is  obtained  by  solution  of  equations  (2.13a-d).  This  system  is  then 
identical  to  that  used  by  Winslow  et  al.'^  and  Carroll  et  al.'^  For  a  number  of  inputs 

P*  (t*)  the  output  from  the  system  is  analytic.  Table  1  below  lists  several  of  these 
solutions  (Carslaw  and  Jaeger^^). 

Table  1 -Sample  PSP  Indicated  Pressures 


Pq  +P,*step(t') 


Po*  +P,*sin(a)'t*) 


exp 


\2  71 


(2i  +  iy-t 
'  4 


.      .  tanh(p)  .  /  .  .\  „     j—r  . 
Po  +Pi  — ^sm(co  t  j,  P  =  Vjco  ,  CO  = 


Note  that  the  function  step(t),  sometimes  referred  to  as  the  unit  step  function,  is  zero  for 
t  <  Oand  unity  for  t  >  0.  Also  note  that  the  solution  shown  for  the  sinusoidal  input  is 
only  that  due  to  the  steady  periodic  conditions,  not  the  initial  transient  period. 

Shown  in  Figure  2.3  is  the  simulated  step  response  using  the  fu-st  1, 2, 4,  and  20 
terms  from  the  analytical  solution  of  Table  1.  Labeled  on  the  graph  are  several  points 
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Figure  2.3--Simulated  PSP  unit  step  response  using  diffusion-based  model  with  linear 
calibration.  Notes:  A)  t;,^  =  1.1,  B)  t;,,  =  I2/71' ,  C)  t;^  =1.8. 
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Figure  2.4--Simulated  PSP  frequency  response  using  diffusion-based  model  with  a  linear 
calibration. 
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Figure  2.5--Comparison  of  PSP  data  and  diffusion-based  model  output  with  a  linear 
calibration. 

throughout  the  response.  It  can  be  seen  that  the  time  to  reach  99%  of  the  total  response  is 
approximately  1.80t  as  compared  to  the  value  of  (1 2/71^)1  put  forth  in  the  paper  by 
Borovoy  et  al.'l  It  can  also  be  seen  that  a  model  order  of  N  =  5  should  be  sufficient  to 
represent  most  of  the  dynamic  behavior. 

Figure  2.4  shows  the  simulated  frequency  response  of  the  PSP.  There  are  two 
interesting  items  of  note  here.  First,  the  amplitude  slopes  off  at  -10  dB  per  decade. 
Second,  the  phase  shift  levels  off  at  ^5°.  Following  the  trend  of  l"  and  2"**  order 
systems  it  would  then  seem  that  this  system  exhibits  the  behavior  of  a  '"V"  order  system, 
which  has  some  meaning  in  the  context  of  fractional  calculus  in  which  derivatives  and 
integrals  of  rational  order  are  defined.  In  the  text  by  Oldham  and  Spanier"  on  this 
subject,  the  relation  between  many  diffiision-based  phenomena  and  the  14  order 
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derivative,  or  semiderivative,  is  shown.  An  attempt  has  been  made  to  express  this  current 
model  using  fractional  calculus,  but  no  clear  results  have  yet  been  achieved. 

Figure  2.5  shows  a  set  of  experimental  input  and  output  pressure  data.  Also 
included  in  the  figure  is  the  model  output  pressure  based  on  the  linearized  diffusion 
model  of  equations  (2.25a,b)  with  N  =  5 .  The  value  of  the  time  constant  that  minimizes 
the  RMSE  between  the  PSP  indicated  pressure  and  the  model  output  pressure  is 
T  =  0.437 .  The  RMSE  of  the  curvefit  is  0.91  kPa  which  is  a  very  good  fit.  The  thickness 
of  this  particular  sample  was  17  ^m.  Using  x  =  a^/D  a  value  of  D  =  6.6  x  10"^  cm^/s  is 
acquired  for  the  mass  difftisivity,  which  compares  well  to  values  fi-om  other  sources.^'^'^^ 
Diffusion  Based  Model  with  Stem-Volmer  Calibration 

The  most  widely  used  calibration  for  PSP  is  based  on  the  Stem-Volmer  relation. 
Taking  equations  1.5  and  1.7  and  letting  Q*  =  Q/P„f  we  get 

•■=f(''')=fr^    ,    ■        .  <2.28a) 

or 

P*=g(r)=Q-+(l-Q')(l/r).  (2.28b) 
In  the  above  equations  the  parameter  Q*  uniquely  determines  the  static  calibration. 
Figure  2.6  illustrates  the  basic  difference  between  the  linear  calibration  of  equafion  (2.18) 
and  the  Stem-Volmer  calibration  of  equation  (2.28).  Note  that  this  figure  shows  the 
intensity  versus  pressure,  as  opposed  to  the  inverse  of  intensity  versus  pressure  which  is 
normally  presented.  Applying  equation  (2.8)  to  equations  (2.28a,b)  and  substittiting  into 
equation  (2.16a)  gives 


31 


Pressure 

Figure  2.6~Comparison  of  calibrations. 

dK^  d'K  2        [dK]'  af(p'(t')) 

5t*"(ax7"K+f(p:(t-))U*J  "    a-    •  ^^-^^^ 

This  is  the  governing  equation  for  the  PSP  dynamics  based  on  the  Stem-Volmer 
calibration.  It  can  be  seen  that  the  middle  term  on  the  RHS  is  highly  nonlinear  and 
contains  a  forcing  component  as  well.  Because  of  this  term,  a  numerical  solution  to  this 
equation  would  require  a  complex  numerical  algorithm  to  be  developed.  However,  rather 
than  perform  a  direct  solution  to  this  equation,  a  different  method  will  be  employed.  It 
will  begin  by  solving  for  the  oxygen  concentration  from  equations  (2.13a-d).  These 
equations  can  be  transformed  to  a  linear  partial  differential  equation  with  forcing  with 
homogeneous  boundary  conditions  and  initial  condition  via 

m(x',t')=n*(x',t*)-P.:(t*).  (2.30) 
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Then,  modal  analysis  can  be  applied  (see  equations  (2.20)  and  (2.21)).  The  resulting 
system  equation  is  identical  to  equation  (2.25a).  The  system  equation  together  with  an 
output  equation  designed  to  give  the  oxygen  concentration  distribution  is  then 


r  : 
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0" 

f  : 
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•) 
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\  '• 
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\  '• 

J 

0  -2sm{X.,)/X; 


•>•(.■)■ 

m 


(2.31a) 


n'(x',t")=(--,cos(>.|X'),  ■• 


a,(f) 


(2.31b) 


This  system  is  linear  and  can  be  easily  solved  at  a  discrete  number  of  points  throughout 
the  layer  to  achieve  an  oxygen  concentration  distribution.  This  oxygen  concentration 
distribution  will  then  be  converted  to  an  intensity  distribution  by  using  the  calibration  in  a 
dynamic  sense, 


(2.32) 


n  (x  ,t  j-Q 

It  may  be  helpful  for  the  reader  to  refer  once  again  to  the  discussion  leading  from 
equation  (2.7)  up  to  equation  (2.1 1)  for  the  validity  of  equation  (2.32).  Then,  the 
intensity  distribution  can  be  numerically  integrated  to  get  the  integrated  intensity  output. 
Making  use  of  the  static  calibration  of  equation  (2.28b)  achieves  the  indicated  output 
pressure. 

Figures  2.7  and  2.8  show  the  simulated  response  of  PSP  to  a  positive  and  a 
negative  step  in  pressure  respectively.  For  the  positive  step,  the  initial  pressure  was 
Pq*  =  0.1 ,  the  final  pressure  was  P*  =  1.0  and  the  value  of  Q*  was  Q*  =  -0.25 .  For  the 
negative  step  the  values  of  Pj  and  P,*  were  reversed.  Shown  in  the  two  figures  are 


comparisons  of  the  two  models  based  on  the  Hnear  caHbration  and  the  Stem-Volmer 
caHbration.  For  the  datasets  based  on  the  hnear  caHbration,  a  model  order  of  N  =  5  was 
used.  For  the  data  based  on  the  Stem-Volmer  calibration,  a  model  order  of  N  =  5  was 
used  to  calculate  the  oxygen  concentration  distribution.  This  was  evaluated  at  21  evenly 
spaced  grid  points  throughout  the  depth  of  the  layer.  Next,  it  was  converted  to  intensity 
per  unit  depth  via  the  calibration  equation  and  numerically  integrated  across  the  layer 
using  the  trapezoidal  rule  integrator.  The  integrated  intensity  was  then  converted  to  an 
indicated  pressure  by  means  of  the  calibration.  From  Figure  2.7  it  can  be  seen  that  for  a 
positive  step  in  pressure,  the  linear  model  solution  leads  the  Stem-Volmer  model 
solution.  In  Figure  2.8  this  behavior  is  reversed.  There  are  several  items  of  note  with 
respect  to  these  comparisons.  First,  the  difference  between  the  models  is  a  function  of 
each  of  the  values  Pj,  P'  and  Q* .  If  the  pressure  data  were  scaled  via  the  ratio 
(P*  -  Pq  )/(Pi*  -  Pq*  )  the  response  of  the  Stem-Volmer  based  model  would  continue  to 
show  different  behavior  for  different  values  of  Pj ,  P,*,  and  Q*  whereas  the  linear 
calibration  based  model  would  always  show  the  same  response.  Second,  as  AP  of  the 
step  decreases,  the  two  models  will  collapse  to  the  same  behavior.  Similarly,  the  models 
will  also  collapse  to  the  same  behavior  for  Q'     -oo . 

Figure  2.9  shows  experimental  input  and  output  data  for  a  positive  step  in 
pressure.  Also  shown  is  a  comparison  of  curvefits  of  the  two  preceding  diffusion-based 
models.  The  results  of  the  curvefits  are  shown  in  Table  2.  For  the  model  based  on  the 
Stem-Volmer  calibration,  a  value  of  Q*  =  -0.163  was  used.  It  is  quite  clear  that  the 
model  using  the  Stem-Volmer  calibration  does  a  much  better  job  of  approximating  the 
PSP  response  than  the  model  based  on  the  linear  calibration. 


Table  2~Comparison  of  Curvefit  Values 

Calibration 

T  (sec) 

RMSE  (kPa) 

Linear 

1.25 

4.10 

Stem-Volmer 

1.11 

1.39 

Figure  2.7--Simulated  positive  step  in  pressure  (model  comparison). 
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Figure  2.8--Simulated  negative  step  in  pressure  (model  comparison). 
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Figure  2.9--Positive  step  in  pressure  from  experimental  data  (model  comparison). 


CHAPTER  3 
DYNAMIC  COMPENSATION 

The  task  of  designing  a  compensation  system  for  PSP  is  comprised  of  two  steps. 
First,  a  model  must  be  developed  for  the  response  of  the  PSP  to  arbitrary  dynamic 
pressures.  Such  models  for  the  PSP  dynamics  have  already  been  presented  in  the 
previous  chapter.  The  second  step  makes  use  of  the  PSP  model  in  combination  with 
feedback  control  to  generate  a  compensated  output.  This  chapter  will  first  show  a 
complete  presentation  of  the  general  compensation  scheme  known  as  Model  Reference 
Control  (MRC).  Then,  it  will  progress  into  the  specifics  of  the  application  of  MRC  to 
PSP.  Finally,  the  chapter  will  wrap  up  with  theoretical  simulations  and  compensation  of 
some  of  the  experimental  data  presented  in  the  previous  chapter. 

Model  Reference  Control 

Model  Reference  Control  is  a  generalized  scheme  in  which  one  can  take  a 
dynamical  system  (plant,  Gp)  with  known  or  modeled  parameters  and  with  the  use  of 
feedback  control,  give  the  resulting  system  (model  reference,  Gm)  a  desired  type  of 
dynamics.  The  system  shown  in  Figure  3.1  is  similar  to  that  presented  in  standard 
Adaptive  Controls  texts  with  regard  to  MRC  with  two  exceptions. First,  there  is  not  a 
constant  gain  feedback  term  from  the  output  Y.  It  is  stated  in  an  exercise  later  in  the 
presentation  of  loannou  and  Sun  that  such  a  term  is  unnecessary.  Second  is  the  addition 
of  the  gain  1/kp  in  the  feedback  from  Y.  This  addition  was  found  to  simplify  some  of  the 
later  analysis. 
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Figure  3.1 --Structure  of  the  Model  Reference  Control  scheme. 

The  MRC  design  begins  with  a  set  of  assumptions  (either  known  or  given):''^ 
Al .      Gp  (s)  =  kpNp  (s)/Dp  (s)  is  a  known  transfer  function. 

A2.  Np(s)  is  a  monic  Hurwitz  polynomial  of  degree  np. 

A3.  Dp(s)  is  a  monic  polynomial  of  degree  dp. 

A4.  The  relative  degree  of  Gp(s)  is      =  dp  -  np  _ 

A5.  Gm  (s)  =  k„,N^  (s)/D„  (s)    a  user-defined  transfer  function. 

A6.     Nm(s)  and  Dm(s)  are  monic  Hurwitz  polynomials  of  degree  Um  and  dm 
respectively,  where      -  dp . 

A7.     The  relative  degree  of  Gni(s)  is      =  d„  -  n^  ^  ^here  ^m=\. 
In  the  above  assumptions,  the  term  "monic"  means  that  the  coefficient  of  the  highest 
order  polynomial  term  is  unity.  The  term  "Hurwitz"  implies  that  the  roots  of  the 
specified  polynomial  all  lie  in  the  lef^  complex  half-plane.  As  a  result  of  the  assumptions 
(Al-7),  we  will  also  define  H(s)  =  N(s)/D(s) ,  where  N(s)  =  (s"'"',  s'''-\ ... ,  l)^  and  D(s) 
is  an  as-yet  undetermined  monic  polynomial  of  degree  d  =  dp .  Notg  that  H(s)  is  a  single- 
input,  multi-output  system.  Analyzing  the  system  of  Figure  3.1  yields 
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°  "   "  (3.1) 


U  DDp-DpC,N-NpC,-N 


where  the  Laplace  variable,  s,  has  been  dropped  from  the  notation  for  simplicity.  It  is 
desired  by  the  user  for  the  Model  Reference  System  (Closed-Loop  System)  to  have  the 
dynamics 


Y    ,  N^(s) 

77  =  ^^«.TrTl-  (3-2) 


U  "Djs) 


Equating  the  systems  gives 


CokpNpD  k„N, 


  m  m 


DDp-D^C,.N-N^C,N  D„ 


(3.3) 


An  inspection  of  the  polynomials  on  the  left-hand  and  right-hand  sides  reveals  that  the 
numerators  and  denominators  are  all  monic  polynomials  excepting  the  values  Co,  kp  and 
km.  Therefore,  equating  these  constant  gains  gives 


Cm 
0  =  —  •  (3.4) 

If  we  now  let 


kp 


D  =  DoN„,  (3.5) 
then  equation  (3.3)  can  be  re-arranged  into  the  form 

DpC,  •  N  +  NpC,  •  N  =  D„ (n^D^  -  N^D„ ).  (3.6) 
This  equation  will  be  termed  the  "Matching  Equation"  and  will  be  used  to  calculate  the 
elements  of  the  unknown  vectors  Ci  and  C2.  Before  proceeding  to  this  solution  we  must 
first  define  the  polynomial  Dq.  In  general  it  is  best  to  choose  the  roots  of  Do  to  be  on  the 


39 


order  of  the  roots  of  Dm  to  keep  the  feedback  gains  at  reasonable  levels.  With  all  of  the 
variables  in  equation  (3.6)  having  now  been  defined  except  for  Ci  and  C2,  we  can 
proceed  to  a  solution.  For  the  analysis  presented  here,  we  shall  assume  that  we  are 
dealing  with  a  plant,  Gp,  with  the  property  np  =  dp  - 1 ,  or,  Ap  =  1 .  Assumption  (A7)  of 

the  MRC  scheme  requires  then  that  =  1 .  The  simplest  dynamical  system  that  has  this 
property  is  a  1^'  order  system 


T  s  +  X 

where  T  is  the  time  constant  of  the  desired  system. 

Let  us  for  a  moment  consider  a  plant  with  dp  =  3 .  Let  us  also  assume  the 
following  notations 

C,  -N  =  a,s^  +a,s  +  a. 


(3.7) 


N  =Y2s'  +Y,s  +  Yo 


Dp  =63S  +52S  +5,s  +  5o 


(3.8a-e) 


Do(N„Dp  -NpD„)=S5s'  +ey+s,s'  +S2S'  +8,s  +  eo 

Ifwe  now  define  x  =  (a2,a,,ao,P2,p,,pJ^  and  b  =  (85,e4,e3,82,e,,eo)'",  the 
Matching  Equation  can  be  transformed  into  a  set  of  standard  linear  equations  ( Ax  =  b ) 


"S3 
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0 

0 

0 

0" 
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(3.9) 


Some  notes  regarding  equation  (3.9): 
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B 1 .     Column  1  of  matrix  A  is  simply  the  vector  of  coefficients  of  the  polynomial  Dp(s) 
with  trailing  zeros  added  to  bring  the  length  of  the  vector  up  to  2dp. 

B2.     Column  2  of  A  is  the  same  as  column  1  but  with  the  last  element  rotated  around  to 
the  top  of  the  vector.  Similariy  for  column  3  with  respect  to  column  2. 

33.     Column  4  of  A  is  the  vector  of  coefficients  of  the  polynomial  Np(s)  with  one 

leading  zero  and  with  trailing  zeros  added  to  bring  the  length  of  the  vector  up  to 
2dp. 

B4.     Columns  5  and  6  of  A  are  1 -element  rotations  of  columns  4  and  5  respectively. 
From  this  example,  the  general  pattern  for  the  matrix  A  and  vectors  x  and  b  for  a  plant  of 
degree  dp  can  be  extrapolated.  Solving  this  system  of  linear  equations  gives  the  values  of 
the  tti's  and  Pi's  which  are  in  turn  the  elements  of  the  vectors  Ci  and  C2. 

Application  of  MRC  to  PSP  Svstem 

The  purpose  for  which  the  Model  Reference  Control  scheme  will  be  used,  with 
regard  to  PSP,  is  not  one  for  which  it  was  intended.  That  is,  the  MRC  scheme  will  be 
used  as  a  "State  Observer".  Considering  a  PSP  plant  with  input  P(t)  and  output  Ppgp  (t), 
we  desire  to  have  a  measurement  of  the  unknown  input  as  opposed  to  what  we  can 
measure  which  is  the  output.  Therefore,  we  have  designed  the  MRC  so  that  the  desired 
measurement  is  an  internal  state.  We  will  further  specify  parameters  of  the  design  so  that 
the  estimated  input  to  the  plant  will  have  very  fast  dynamics  so  that  it  should  quickly 
approach  the  true  pressure  P(t). 

For  the  particular  application  of  PSP  dynamics,  the  MRC  scheme  shown  in  Figure 
3. 1  will  be  changed  to  the  diagram  shown  in  Figure  3.2.  We  will  consider  the  measured 
PSP  pressure,  Pp^p  (t),  to  be  the  input  U.  The  internal  state  V  will  be  the  estimation  of  the 
true  pressure,  P(t)  =  P^  (t) ,  and  will  in  turn  be  our  value  of  the  compensated  pressure. 
The  output  Y  of  the  MRC  scheme  is  simply  an  estimation,  Pp^p  (t),  of  the  input.  The 
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Ppsp 


Figure  3.2--Structure  of  the  Model  Reference  Control  scheme  for  use  as  a 
compensator  for  PSP  data. 

other  internal  states  Wi,  W2,  and  Z  are  of  no  relevance.  The  plant,  Gp,  is  one  of  the 
dynamics  models  for  the  PSP  given  in  Chapter  2.  Our  aim  is  to  design  the  MRC 

parameters  so  that  Gm  is  a  very  fast  system.  Thus,  the  output  Ppsp  (0  will  be  closely 

mimicking  the  input  Ppsp  W.  Furthermore,  since  Ppsp(t)  is  the  output  of  the  plant  Gp,  the 

input  to  the  plant  should  be  closely  tracking  the  true  pressure.  That  is,  ^  =  H^)  =  Pc  0). 
So,  by  picking  off  the  input  to  the  plant  Gp  and  making  it  an  output,  we  will  be  able  to  get 
an  approximate  measure  of  the  true  pressure  versus  time. 

To  begin,  we  will  need  to  use  a  linearized  model  for  the  PSP  dynamics.  It  is 
necessary  to  use  a  linearized  model  in  order  to  calculate  the  values  Co,  Ci,  and  C2.  After 
calculating  the  parameters  of  the  MRC  using  the  linearized  plant,  the  nonlinear  plant  will 
be  resubstituted  into  the  overall  scheme  keeping  the  same  values  for  Co,  Ci,  and  C2.  For 
the  linearized  model  we  will  use  the  diffusion  model  from  Chapter  2  with  the  linear 
calibration  (equations  2.25a,b).  Further,  for  the  purposes  of  presentation,  we  will  assume 


a  2"**  order  model.  The  linearized  plant  specifications  are 
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P(t) 

P(.) 
P(t) 

P(t)J 

Converting  this  from  State-Space  form  into  Transfer  Function  form  we  get 

's'  +  24.67s +  54.79' 


"ct,' 

■-2.47  0.00" 

"a," 

"0.00-1.27" 

+ 

0.00-22.21 

0.00  0.42 

*PSP  0 

)  = 

[  0.64  -0.21] 

"a," 

+ 

[1.00  0.00] 

«2. 

(3.10a,b) 


G(s)  = 


-0.90s -18.22 


(3.11) 


+  24.67s +  54.79 

The  numerator  of  this  equation  is  an  array  with  two  rows  -  one  for  each  of  the  inputs 
P(t)  and  P(t).  Notice  that  the  transfer  function  for  P(t)  is  that  of  an  all-pass  system. 
Making  use  of  the  fact  that  P(t)  =  s  •  P(t),  we  can  reduce  the  system  of  equation  (3.1 1) 
from  a  MISO  system  (multi-input,  single-output)  to  a  SISO  system  (single-input,  single- 
output) 

G,(s)=0.0994(;^\Q-%Xs^54.88L  (  10.05,-54.88) 

(s  +  2.47Xs  + 22.21)  (-2.47,-22.21)  ^  ^ 

The  plant  expressed  in  this  form  has  an  input  of  P(t)  and  output  of  Ppgp  (t).  This 
equation  will  be  used  to  define  the  plant  during  the  design  stage  of  the  MRC 
development.  Later,  the  nonlinear  form  of  the  plant  (diffiision  model  with  Stem-Volmer 
calibration)  will  be  used  in  its  place.  The  latter  portion  of  equation  (3.12)  includes  a 
definition  of  the  plant  in  "Zero-Pole"  format  where  the  zeros  and  poles  of  the  transfer 
function  are  listed  in  vector  form.  The  value  of  0.0994  in  front  is  the  high-frequency  gain 
of  the  fransfer  function. 

In  the  earlier  section  of  this  chapter,  the  MRC  scheme  was  developed  for  systems 
with  the  property  Ap  =  1 .  Notice  however,  that  the  plant  definition  of  equation  (3.12)  has 
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a  relative  order  of  Ap  =  0 .  Therefore,  we  must  either  re-develop  the  MRC  scheme  for 
systems  with  Ap  =  0 ,  or  find  a  way  to  modify  the  plant  to  give  it  a  relative  order  of 
Ap  =  1 .  There  is  a  simple  "trick"  that  can  be  employed  to  give  the  desired  effect  of 
changing  the  relative  order  of  the  plant.  By  augmenting  the  plant  with  a  very  fast  1^' 
order  system,  the  resulting  transfer  function  will  have  the  desired  properties.  It  should  be 
noted  that  this  augmented  1^'  order  system  will  also  have  to  operate  on  the  actual  PSP 
data  prior  to  becoming  the  input  of  the  MRC  scheme.  The  question  now  is  how  to  choose 
the  time  constant  of  the  augmented  system.  Inspection  of  the  poles  of  the  plant  shows 


poles  =  - 


(2i-l)- 
'l 


2 


,i  =  l,2,...  (3.13) 


Following  this  pattern,  the  next  value  in  the  sequence  will  be  used  as  the  pole  of  the 
augmented  system  yielding: 

G  (s)=k  ^^i^  k  (s +  10.05X5  + 54.88) 

In  using  this  choice  of  augmentation  there  will  be  no  loss  of  dynamics  fi-om  the  physical 
system.  That  is,  because  the  augmented  system  is  in  essence  a  lowpass  filter,  but  with  a 
very  high  cutoff  fi-equency,  the  dynamics  of  the  physical  system  will  be  relatively 
unaffected.  Considering  this  augmented  plant  as  the  design  plant  (Gp  =  Gpa)  for  the 
purposes  of  calculating  the  Model  Reference  Controller  parameters,  we  can  choose  the 
model  system  as 

Gjs)=k,^  =  100.^  (3.15) 
D^(s)  s  +  100  ' 
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which  is  a  l"  order  system  with  a  time  constant  of  0.01.  For  the  Plant  Gpa(s)  of  equation 
(3.14)  and  Model  Gm(s)  of  equation  (3.15)  and  for  the  choice  of  Do(s) 

Do(s)  =  (s  +  lOOXs  + 1  lOXs  + 120),  (3.16) 
the  solution  to  the  Matching  Equation  gives 

Co  =16.32 

C,  =(-7.86el,  1.33e4,  1.13e6)  (3.17a-c) 
=(-3.79e4,  -4.39e6,  -1.3 leg) 

Or,  if  we  redefine  the  feedback  systems  in  the  following  manner: 


H(3)-C.-N(s)  N,(s) 

H,(s)=^  =  k.^ 
D(s)  'D,(s) 


(3.18a,b) 


we  get  (in  zero-pole  format) 


H,(s)= -78.57  ,t^?:2MMlL 
(-100,-110,-120) 

„,(s).-37,47Jr^Z™L 
(-100,-110,-120) 

Following  this  procedure,  the  MRC  parameters  can  be  calculated  for  a  PSP  plant  of 
arbitrary  order. 

MRC  Simulations  and  Compensation  of  Experimental  Data 
It  now  remains  to  be  seen  how  well  the  MRC  scheme  can  be  used  to  compensate 
the  PSP  data.  Using  Matlab  and  Simulink,  programs  were  developed  to  implement  the 
MRC  scheme  (see  Appendices  A,  B).  With  this  code  several  simulations  were 
performed.  In  these  simulations,  a  5'^  order  plant  design  and  its  accompanying  MRC  was 
used  as  opposed  to  the  2"''  order  design  already  presented.  There  were  three  basic  sets  of 
tests  performed.  The  first  used  the  linear  form  of  the  Plant  to  create  the  input  Ppsp(t),  and 
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Figure  3.3~Linear  compensation  simulation  of  a  positive  step  in  pressure. 


the  linear  form  of  the  Plant  as  Gp(s).  This  test  was  performed  to  ensure  that  the 
compensator  works  as  it  was  designed,  that  is,  the  compensated  output  should  have  the 
dynamics  of  the  ideal  design  Model  Reference  System.  The  performance  of  this  linear 
version  of  the  MRC  compensator  was  tested  on  a  positive  step  from  0.1  atm  up  to  l.Oatm. 
The  results  of  this  test  can  be  seen  in  Figure  3.3.  Notice  that  the  PSP  indicated  pressure 
takes  approximately  two  time  units  to  come  to  equilibrium,  whereas  the  compensated 
output  levels  out  in  less  than  0.05  time  units.  Therefore,  this  form  of  compensation 
system  does  indeed  work  in  the  linear  case  and  affects  a  more  than  forty-fold  increase  in 
speed.  Thus,  the  MRC  compensation  scheme  does  work  in  the  linear  limit.  The  second 
set  of  tests  used  the  full  nonlinear  form  (diffusion  model  with  Stem-Volmer  calibration) 
to  create  the  input  and  as  Gp(s).  The  compensator  parameters,  Co,  €,,  and  C2,  used  in  this 
test  were  from  the  linear  design.  The  response  of  the  compensator  system  to  a  positive 
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Figure  3.4--Compensation  simulations  using  nonlinear  from  of  the  PSP  plant,  (a)  Positive 
step  in  pressure;  (b)  Negative  step  in  pressure. 
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Figure  3.5--Compensation  simulations  using  nonlinear  form  of  the  PSP  plant  with  a 
sinusoidal  driving  pressure,  (a)  to*  =  1 ;  (b)  co*  =  3 ;  (c)  ©*  =  10 ;  (d)  to*  =  30 . 
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Figure  3.6--Compensation  of  experimental  data,  (a)  Positive  step  in  pressure;  (b) 
Sawtooth  driving  pressure. 
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and  negative  step  was  simulated.  As  in  Chapter  2,  the  step  simulations  were  between  the 
pressures  0.1  atm  and  l.Oatm  and  Q  =  -0.25  determined  the  Stem-Volmer  calibration. 
Shown  in  Figure  3.4  are  the  results  of  the  simulations  for  the  positive  step  and  negative 
step.  The  compensated  response  to  a  positive  step  in  pressure  contains  sizable 
oscillations,  whereas  the  negative  step  response  is  practically  that  of  the  ideal  linear 
compensator.  Thus,  the  nonlinearity  of  the  PSP  system,  as  it  pertains  to  the  MRC 
compensator,  is  more  in  evidence  for  positive  changes  in  pressure  than  for  negative 
changes.  The  third  set  of  simulations  were  similar  to  the  second  in  that  the  nonlinear 
form  was  used  for  the  plant.  However,  rather  than  applying  impulsive  steps  in  pressure, 
sinusoidal  fluctuations  in  pressure  were  simulated.  The  sinusoidal  tests  were  between  the 
limits  0.1  atm  and  l.Oatm  with  the  same  value  of  Q.  The  frequency  of  the  driving 
pressure  was  set  at  four  different  values,  to*  =  1, 3, 10,  and  30 .  The  results  of  these  tests 
are  shown  in  Figure  3.5.  There  are  several  items  to  note  concerning  these  simulations.  It 
can  be  seen  that  the  compensator  does  a  remarkable  job  even  out  to  driving  frequencies 
ten  times  the  critical  frequency  of  the  paint.  One  thing  of  note  with  respect  to  the 
sinusoidal  simulations  is  that  the  phase  shift  of  the  compensated  data  relative  to  the 
applied  pressure  is  very  nearly  0.01  in  each  of  the  four  cases  (frequencies)  shown.  Thus, 
the  compensator  is  nearly  a  "linear  phase"  system.  Also  note  that  this  value  of  0.01  is 
exactly  the  time  constant  of  the  designed  Model  Reference  System. 

What  remains  is  to  determine  how  well  the  compensator  works  on  actual  PSP 
data.  Figure  3.6  shows  the  compensator  response  from  two  different  experiments.  The 
first  is  a  positive  step  in  pressure  from  20  kPa  up  to  120  kPa  (note:  Pref  =  100  kPa).  The 
time  constant  of  the  PSP  was  1.1 1  sec.  The  second  was  a  sawtooth  driving  pressure  and 


the  time  constant  was  0.437  sec.  In  the  step  test  a  shght  overshoot  can  be  seen,  but  the 
compensator  is  still  a  great  improvement  over  the  raw  PSP  data.  The  compensation  of 
the  sawtooth  data  has  a  peak  error  of  4  kPa  and  an  RMS  error  of  about  1  kPa.  The 
maximum  error  occurs  on  the  rising  edge  of  the  sawtooth  wave  where  the  highest 
frequencies  are  contained.  From  these  two  sets  of  experimental  data,  we  can  conclude 
that  the  MRC  compensation  scheme  does  indeed  work.  The  next  step  will  be  to  apply 
this  to  a  sequence  of  full-field  images. 


CHAPTER  4 
EXPERIMENTAL  SETUP  AND  PROCEDURES 

This  chapter  will  begin  with  a  presentation  of  the  two  facilities  employed  to 

acquire  the  experimental  data  presented  in  this  dissertation.  Both  of  these  facilities  are 

located  at  the  University  of  Florida,  Department  of  Aerospace  Engineering,  Mechanics 

and  Engineering  Sciences.  The  first  of  these  is  a  test  cell  made  for  testing  the  dynamic 

response  of  small-sized  PSP  samples.  The  second  is  a  test  stand  for  investigation  of  a 

pulsed  transverse  jet  in  a  Mach  1.56  crossflow.  The  end  of  this  chapter  will  be  a 

presentation  of  the  image  processing  techniques  used  to  get  a  single  calibrated  PSP 

image.  The  method  of  compensating  these  calibrated  images  will  be  presented  in  the 

following  chapter. 

PSP  Dynamics  Test  Cell 
The  data  presented  in  Chapters  2  and  3  were  performed  in  the  PSP  Dynamic  Test 
Cell.  This  test  cell  is  designed  to  control  the  pressure  of  the  air  over  the  PSP,  the 
temperature  of  the  PSP  sample,  and  the  temperature  of  the  air  entering  the  test  cell.  The 
experimental  apparatus  consists  of  the  test  cell  itself  (with  various  accessories),  an 
external  heat  exchanger,  a  lamp,  a  photomultiplier  tube  (PMT),  some  optical  components, 
a  computer  with  data  acquisition  card  and  various  intermediate  circuits.  Additionally, 
vacuum  and  high  pressure  lines  were  available  to  connect  to  the  test  section.  The  entire 
apparatus  rested  on  an  optical  bench  with  an  enclosing  hood  to  shut  out  external  light 
sources. 
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A  front  view  of  the  test  cell  is  shown  in  Figure  4.1.  A  schematic  view  of  the  test 
cell  is  also  shown  in  Figure  4.2.  There  are  six  ports  (1/8"  NPT)  around  the  periphery  of 
the  cell  to  allow  access  to  the  small  volume  (1/8"  thick  x  1/2"  diameter)  immediately 
above  the  sample.  These  ports  can  be  instrumented  in  a  variety  of  ways.  As  shown  in  the 
Figure  4. 1,  it  is  instrumented  with  two  digital  valves  (Asco/Angar  AL41 12L),  a  high 
frequency  response  pressure  transducer  (Entran  EPV-501X-25AZ),  and  a  thermistor 
(Omega  44201).  The  two  remaining  ports  are  capped.  The  two  digital  valves  are  hard- 
wired to  be  180°  out  of  phase  so  the  test  cell  is  either  open  to  vacuum  or  open  to  high 
pressure.  With  this  feature,  a  pressure  step  from  vacuum  to  ambient  pressure  occurs  in 
just  under  20ms.  Additionally,  an  analog  valve  (Pneufronics  VSO,  1.25mm  orifice,  not 
shown  in  figure)  can  be  mounted  to  one  of  the  ports.  With  the  use  of  the  analog  valve, 
the  pressure  in  the  cell  can  be  driven  at  any  arbitrary  function  the  user  may  desire  such  as 
a  sine  wave,  square  wave  or  sawtooth  wave.  It  should  be  noted  that  because  of  the  small 
orifice  of  the  analog  valve,  the  pressure  inside  the  cell  was  band-limited  to  less  than 
80Hz.  However,  with  the  dual  digital  valves  confrolling  the  pressure  in  the  cell, 
frequency  content  upwards  of  250Hz  could  be  achieved,  but  only  in  the  form  of  a  step  or 
square  wave.  It  was  later  discovered  that  using  a  latex  burst  diaphragm  on  one  of  the 
ports  could  create  frequency  content  of  the  pressure  signal  in  excess  of  IkHz,  however, 
no  such  experiments  are  presented  in  this  dissertation.  The  window  on  the  front  of  the 
cell  is  plexiglass  and  allows  optical  access  to  the  PS?  sample.  The  sample  is  a  small 
piece  of  aluminum  (12.7mm  x  15.9mm  x  3.2mm)  with  the  PSP  applied  to  one  of  the 
faces.  The  size  of  the  samples  was  chosen  so  they  could  be  easily  and  cheaply  made 
from  12.7mm  x  3.2mm  (1/2"  x  1/8")  aluminum  bar  stock.  The  sample  is  mounted  on  top 
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of  a  raised  circular  platform.  An  0-ring  around  the  base  of  the  platform  seals  the  interior 
of  the  cell.  Behind  the  platform  is  a  double  stack  of  thermoelectric  chillers  (TEC's) 
(Melcor  CP1.0-31-05L)  to  heat  or  cool  the  sample.  Mounted  further  behind  the  TEC's 
are  a  set  of  cooling  fins  and  a  small  DC  fan.  Originally  the  test  cell  was  made  with  a 
water-bath  behind  the  circular  platform  to  control  the  temperature  of  the  sample,  but  this 
was  later  changed  out  for  the  TEC's. 

The  inflow  air  heat  exchanger  unit,  shown  in  Figure  4.3,  is  made  of  two 
aluminum  plates  (17.8cm  x  15.2cm)  clamped  together.  Machined  into  one  of  the  plates  is 
a  zig-zagging  channel  thru  which  air  can  flow.  Mounted  on  top  of  the  the  unit  are  four 
TEC's  (Melcor  CP  1.4-1 27- lOL)  with  cooling  fins  and  fans  to  control  the  temperature  of 
the  top  plate.  The  entire  unit  is  encased  in  foam  insulation  to  aid  in  temperature  control. 
Mounted  at  the  exit  is  a  thermistor  (Omega  44201)  to  measure  the  temperature  of  the  air 
leaving  the  unit. 

Two  different  styles  of  lamp  can  be  used  to  illuminate  the  PSP  sample.  The  first, 
shown  in  Figure  4.4,  is  a  quartz-tungsten-halogen  lamp  with  stabilized  power  supply  and 
photofeedback  system  (Oriel  QTH  lamp  housing  with  68830,  68850,  68855  accessories). 
The  second  style  of  lamp  is  shown  in  Figure  4.5  and  incorporates  an  array  of  super-bright 
LED's  (ISSI LM2). 

The  PMT  used  to  collect  and  measure  the  light  fi-om  the  PSP  sample  has  a 
spectral  response  from  185nm  to  900nm  and  a  bandwidth  of  20kHz  (Hamamatsu  HC120- 
05),  although  analog  RC  lowpass  filters  used  to  aid  in  reduction  of  shot  noise  lessened  the 
bandwidth. 


There  are  several  components  of  optics  in  the  setup  -  two  filters,  a  mirror  and  a 
lens.  The  first  of  the  two  filters  is  a  bandpass  filter  centered  at  450nm  (Melles  Griot 
03FIB004)  and  is  mounted  to  the  lamp.  The  second  filter  is  also  a  bandpass  filter,  but  is 
centered  at  650nm  (Melles  Griot  03FIB014)  and  is  mounted  to  the  fi-ont  of  the  PMT.  The 
mirror  reflects  the  excitation  fi-om  the  lamp  into  the  test  cell  and  onto  the  PSP  sample. 
The  light  emitted  from  the  sample  is  then  focused  thru  a  lens  into  the  PMT.  A  schematic 
for  the  path  of  the  light  is  shown  in  Figure  4.6. 

The  computer  used  in  the  experiment  was  one  of  several  IBM-compatible 
Pentium  or  Pentium  II  machines.  A  data  acquisition  card  (DAQ  card)  (National 
Instruments  AT-MIO-16E-2  with  SCB-68  accessory)  was  used  to  sample  the  signals 
fi-om  the  various  transducers  and  to  control  the  operation  of  the  digital  and  analog  valves. 
Data  acquisition  and  control  programs  were  developed  using  National  Instruments' 
LabView  programming  language.  This  program  was  designed  in  a  very  "hands  off 
fashion  and  allowed  the  user  to  control  the  digital  and  analog  valves  and  to  perform  the 
data  acquisition  of  the  transducers  all  from  a  seat  in  front  of  the  computer.  Data  from  the 
fransducers  was  acquired  at  rates  on  the  order  of  5kHz  and  then  downsampled  with 
averaging  to  a  rate  of  IkHz. 

Several  circuits  were  used  to  interface  the  fransducers  and  valves  to  the  DAQ 
card.  First  order  analog  lowpass  filters  with  a  cutoff  frequency  of  480Hz  were  applied  to 
the  signals  from  the  pressure  fransducer  and  from  the  PMT  to  aid  in  noise  rejection. 
Because  the  DAQ  card  could  only  source  5mA  of  current,  power  transistors  were  used  to 
provide  current  amplification  to  the  voltage  signals  coming  from  the  analog  output  and 
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digital  output  lines  of  the  DAQ  card  leading  to  the  valves.  Finally,  resistor  circuits  were 
used  to  create  a  linear  voltage  versus  temperature  output  from  the  thermistors. 

Supersonic  Wind  Tunnel 
The  second  facility  used  for  experimentation  was  a  blow-down  wind  tunnel.  Two 
large  tanks  external  to  the  building  provide  an  initial  source  of  approximately  30m^  of  air 
at  a  maximum  pressure  of  1375kPa.  Supplying  air  to  the  tanks  is  a  300hp  rotary  screw 
compressor  (Quincy  QSI-1000)  capable  of  delivering  34kg/min  at  1375kPa  (950SCFM  at 
200psig).  Attached  to  the  compressor  are  dual  dessicant  dryers  (Hydronix)  that  cycle  in 
turn  to  remove  moisture  from  the  air  after  compression.  A  3"  line  runs  from  the  tanks, 
into  laboratory,  and  thru  one  of  two  control  valves  (Fisher  ET-l'/z" ,  Fisher  ET-1"  with 
Type  667  Actuators,  Series  3582  Valve  Positioners  and  Type  67AF  Filter  Regulators) 
depending  on  the  desired  flowrates  or  stagnation  pressures  required  at  the  test  stand. 
These  valves  are  shown  in  Figure  4.7.  The  two  valves  are  mounted  in  parallel  on  the  3" 
line.  The  valves  are  ultimately  user-controlled  in  one  of  two  fashions.  First,  they  may  be 
controlled  manually  via  a  pneumatic  controller  system  (Fisher  "Wizard"  Controller). 
Alternatively,  the  valves  may  be  confrolled  via  a  computer  connected  to  an  external 
transducer  (Fisher  Type  846  Electro-Pneumatic  Transducer).  The  second  of  these  was 
the  method  of  choice.  The  final  element  of  the  general  test  stand  was  a  large  cast-iron 
pressure  plenum,  or,  stagnation  chamber  shown  in  Figure  4.8.  The  air  flow  is  from  left  to 
right  in  the  image.  To  this  chamber  can  be  mounted  whatever  test  section  the  user  may 
design  such  as  the  Mach  1.56  nozzle  used  in  the  experiments  to  be  presented  in  the  next 
chapter. 
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The  test  section  chosen  for  experiments  reported  in  this  dissertation  has  a  design 
Mach  number  of  1 .60  and  its  geometry  was  calculated  using  the  inviscid  nozzle  design 
program  of  Carroll  et  al.'"  Later  pressure  measurements  made  in  the  tunnel  yield  a  Mach 
number  of  1 .56.  This  difference  is  due  to  the  growth  of  the  boundary  layer  downstream 
of  the  sonic  throat.  A  slight  divergence  of  the  test  section  downstream  of  the  nozzle  exit 
could  have  compensated  for  the  boundary  layer  growth,  however,  it  was  decided  to  use  a 
straight  section  to  ensure  a  perpendicular  viewing  of  the  jet  to  be  described  later.  The 
main  section  of  the  tunnel  downstream  of  the  nozzle  is  a  5.08cm  x  5.08cm  (2"  x  2") 
cross-section  with  optical  access  via  a  15.24cm  (6")  long  plexiglass  wall.  There  are  also 
two  opposing  sidewalls  with  interchangeable  plates  with  provisions  for  additional  optical 
access,  or  transducer  measurements,  or  for  placing  obstructions  into  the  flow.  An 
external  view  of  the  test  section  is  shown  in  Figure  4.9.  The  direction  of  airflow  is  from 
left  to  right.  The  main  plexiglass  window  is  toward  the  front.  Along  the  upper 
removable  sidewall  are  an  array  of  static  pressure  ports  and  thermocouples.  For  a  lower 
viewing  angle,  a  4.44cm  (1.75")  diameter  circular  plug  opposite  the  main  plexiglass  wall 
would  also  be  seen.  Figure  4. 10  shows  a  closeup  view  of  the  front  of  the  plug.  For  the 
purposes  of  the  experiments  to  be  presented  in  the  next  chapter,  this  plug  has  been 
modified  to  create  a  fransverse,  or  normal,  jet  3mm  (0. 1 1 8")  inches  in  diameter 
surrounded  by  four  0.5mm  diameter  (0.020"  0)  pressure  taps  at  90°  increments  on  a 
4.5mm  (0. 1 77")  radius  about  the  jet.  Additionally,  a  single  thermocouple  is  imbedded  ' 
2.5mm  (0. 1 ")  from  the  interior  surface  of  the  test  section.  Figure  4. 1 1  shows  a  rear  view 
of  the  plug.  The  radial  location  of  the  thermocouple  can  be  seen.  Also  shown  in  the 
figure  is  the  digital  valve  used  to  control  the  jet.  The  diameter  of  the  jet  was  chosen  to  be 
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approximately  equal  to  the  free-stream  boundary  layer  thickness  as  calculated  by 
Carroll.''^  Other  studies  have  used  this  same  method  to  choose  the  diameter  of  a  jet 
injected  normally  into  a  supersonic  air  stream.^'^'  The  internal  geometry  of  the  jet  is 
simply  a  straight  3mm  diameter  cylinder.  Other  studies  have  used  a  converging  geometry 
at  the  jet  orifice  to  ensure  a  sonic  throat.  However,  because  the  valve  is  being  cycled,  and 
because  it  is  desired  to  have  as  fast  responding  a  jet  as  possible,  the  minimal  volume 
straight  cylinder  geometry  was  chosen. 

The  optical  system  used  to  measure  the  intensity  of  the  PSP  consisted  of  a 
camera,  a  strobe  and  appropriate  filters.  The  camera  was  a  Photometries  Model  CH250 
CCD  camera  as  shown  in  Figure  4.12.  The  CCD  array  itself  was  a  512  x  512  array  of 
pixels.  The  full-well  capacity  of  the  pixels  is  approximately  320,000  electrons  and  the 
pixels  were  read  out  using  a  14-bit  A/D  converter.  The  strobe,  also  shown  in  Figure  4.12, 
was  a  Xenon  DT-301  Digital  Stroboscope  Tachometer.  The  strobe  could  operate  at 
frequencies  of  2Hz  -  200Hz.  The  duration  of  a  single  flash  of  the  strobe  is  on  the  order  of 
30^s.  A  composite  of  blue  Schott  glass  filters  (BG12  -  3mm,  BG39  -  3mm)  was  placed 
on  the  output  of  the  strobe  to  dampen  its  output  at  red  wavelengths  and  yet  allow  the  blue 
excitation  through.  A  Melles  Griot  03FIB014  bandpass  filter  centered  at  650nm  was 
placed  on  the  input  to  the  camera  to  filter  out  the  source  illumination  from  the  strobe. 

Various  conventional  transducers  were  used  in  these  experiments.  The  stagnation 
pressure  was  measured  with  an  Omega  PX302-200AV  pressure  transducer.  This 
measurement  was  used  as  the  process  variable  of  a  PID  algorithm  to  control  the  Fisher 
1 V"  valve  and  for  calculating  the  Mach  number  of  the  flow.  All  of  the  temperature 
measurements  were  made  using  Omega  T-type  thermocouples.  The  cold-junction 


reference  temperature  will  be  described  later.  Static  pressure  measurements  were  made 
with  a  combination  of  two  different  ranges  of  Kulite  piezoresistive  pressure  transducers. 
Two  Kulite  Model  XT- 190-25 A  (25psia  range)  and  two  Kulite  Model  XT-190-50A 
(50psia  range)  transducers  were  available  for  measuring  the  pressures  at  the  various  static 
pressure  taps.  One  of  the  the  50psia  transducers  was  positioned  at  the  sidewall  pressure 
port  closest  to  the  nozzle  exit  and  measured  the  static  pressure  for  purposes  of  calculating 
the  Mach  number.  It  was  located  far  enough  upstream  of  the  injected  jet  to  be  outside  of 
its  influence.  A  second  SOpsia  transducer  was  placed  at  the  pressure  tap  spanwise  from 
the  jet.  The  two  remaining  transducers  with  ranges  of  25psia  were  placed  at  the  taps 
upstream  and  downstream  of  the  jet  orifice.  Thus,  the  two  25 's  were  placed  at  the 
locations  of  highest  interest  and  the  two  50's  were  placed  at  the  less  critical  locations. 
The  Kulite  transducers  mounted  to  the  pressure  taps  on  the  plug  insert  were  offset  from 
the  interior  of  the  test  section  by  a  5.08cm  x  1 .6mm0  (2.0"  x  0.063 "0)  volume.  The 
frequency  response  of  the  pressure  transducer/tap  volume  combination  was  not 
calculated,  but  was  assumed  to  have  a  bandwidth  greater  than  IkHz  which  was  the 
maximum  sampling  rate  used  in  the  experiments  presented.  The  pressure  transducers 
were  calibrated  using  a  NIST-traceable  pressure  standard  (Druck  DPI- 145  Multiftinction 
Pressure  Indicator).  Additionally,  prior  to  each  run,  the  transducers  were  "zero- 
calibrated"  to  match  the  ambient  pressure  indicated  by  the  pressure  standard. 

The  jet  could  be  turned  on  or  off  or  modulated  at  a  given  frequency  via  a  Snap- 
Tite  digital  valve  (2W13W-1NB-A8E6,  24V,  9W,  1/8  orifice).  The  valve  was  capable  of 
operating  at  frequencies  upwards  of  lOOHz  although  at  the  upper  range  the  valve  had  a 
tendency  to  stick  in  the  open  position.  At  frequencies  below  80Hz  a  duty  cycle  of  40% 


on  and  60%  off  seemed  to  work  the  best  in  terms  of  the  pressures  seen  while  running  the 
tunnel  and  on  the  basis  of  "non-stick".  The  pressure  upstream  of  this  valve  was  set  using 
a  regulator  coming  off  the  lOOpsi  house  air  supply  within  the  building.  This  pressure  was 
further  measured  using  the  Druck  pressure  standard. 

The  experiment  was  controlled  via  a  data  acquisition  and  control  program 
developed  using  LabView  and  running  on  a  Pentium  III  computer.  This  program 
controlled  two  separate  DAQ  cards  within  the  computer  (AT-MIO-16E-2;  PCI-6024E). 
The  first  board  measured  the  5  pressure  transducers,  the  wall  thermocouple,  and  the 
voltage  signals  being  applied  to  the  jet  valve  and  the  strobe.  Prior  to  performing  the 
experiment,  a  switch  was  set  on  the  external  connector  block  (SCB-68)  that  connected 
channel  0  of  the  board  to  a  cold-junction  temperature  sensor  on  the  connector  block.  A 
measure  of  the  cold-junction  temperature  was  then  recorded  and  used  later  to  properly 
convert  the  thermocouple  voltage  readings  into  temperature.  After  measuring  the  cold- 
junction,  the  switch  was  reset  so  that  channel  0  could  be  used  to  measure  one  of  the 
external  transducers.  One  analog  output  on  the  first  board  was  used  to  control  the  main 
valve  responsible  for  running  the  tunnel.  Using  PID  control,  the  stagnation  pressure  was 
monitored  and  set  to  a  desired  pressure  of  158kPa  (23psia)  which  was  just  above  the 
pressure  necessary  to  fiilly  start  the  tunnel.  The  second  DAQ  card  was  used  solely  for  its 
two  analog  output  lines.  These  were  used  to  drive  the  jet  valve  and  the  strobe.  The  phase 
of  the  strobe  relative  to  the  jet  could  be  controlled  down  to  1/50""  of  a  cycle  at  a 
frequency  of  lOOHz.  This  allowed  the  strobe  to  be  flashed  consistently  at  a  desired  phase 
relative  to  the  cycling  of  the  jet.  Hence,  the  acquired  PSP  images  were  actually  ensemble 
averages  of  hundreds  to  thousands  of  individual  strobe  flashes. 
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The  final  element  of  the  tests  performed  in  this  test  section  are  the  PSP  and  TSP 
coatings.  The  PSP  consisted  of  the  luminophore  PtTFPP  [Pt(II)  meso-tetrakis(penta- 
fluorophenyl)porphine]  interspersed  in  the  polymer  matrix  VPDMS  [vinyl  polydimethyl- 
siloxane].  The  TSP  was  a  ruthenium  derivative  in  clear-coat  acrylic.  These  formulations 
was  made  and  applied  to  the  plug  insert  by  a  member  of  Dr.  Schanze's  group  at  the 
University  of  Florida  Department  of  Chemistry.  The  coating  thickness  varied  from 
painting  to  painting  but  was  nominally  on  the  range  from  10|im  -  20(im. 

Image  Processing 

Achieving  a  calibrated  PSP  image  is  a  very  complex  process.  Numerous  papers 
have  been  written  on  just  this  subject.  The  general  process  is  shown  in  the  flowchart  of 
Figure  4. 13.  To  begin,  one  needs  two  images  -  one  at  the  unknown  run  condition  and  a 
second  at  known  reference  conditions.  Consider  the  two  image  intensities  to  be  !„,„  and 
Iref  respectively.  The  next  step  is  to  subfract  the  dark  image,  Idrk,  from  each  of  the  run  and 
reference  images.  The  dark  image  is  an  image  of  the  same  duration  as  the  run  and 
reference  images  but  with  all  excitation  sources  turned  off  Next  is  a  flat  field  correction 
which  compensates  for  the  pixel-to-pixel  variation  in  sensitivity  of  the  CCD  array.  This 
correction  takes  the  form  of  a  matrix,  Cfr,  the  size  of  the  CCD  array  with  values  all  near 
1.0.  When  this  matrix  is  multiplied  by  the  given  image,  the  differences  in  the  pixel  gains 
is  cancelled  out.  The  task  of  generating  the  matrix  Cfr  is  an  involved  process  in  and  of 
itself  In  the  process  of  running  the  tunnel,  the  model  of  interest  may  move  slightly 
relative  to  its  position  at  the  reference  image.  In  a  later  step,  the  data  contained  in  the  two 
images  will  need  to  be  ratioed.  Therefore,  it  is  of  extreme  importance  that  the  two 
images  be  perfectly  aligned  prior  to  the  ratioing  step.  Often,  markers  are  placed  on  a 
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model  to  be  used  as  registration  points.  Alternatively,  certain  parts  of  a  model  may  be 
used  as  registration  points.  For  example,  oftentimes  a  model  will  have  pressure  taps  in 
addition  to  the  PSP  to  be  able  to  get  very  acurate  measurements  at  a  distrete  number  of 
points.  In  a  PSP  image  these  taps  will  show  up  as  tiny  dots  of  low  intensity.  If  a  number 
of  these  are  scattered  over  the  model,  they  can  be  used  to  calculate  the  movement  of  the 
model  relative  to  its  reference  position.  For  certain  models,  a  simple  pixel-shifting  is 
sufficient.  However,  if  the  model  rotates  relative  to  its  reference  or  warps  due  to  the 
loads  it  experiences,  the  registration  algorithm  must  become  much  more  complex.  For 
the  model  used  in  this  dissertation,  warp  and  rotation  are  not  present.  So,  a  simple  linear 
algorithm  will  be  used  to  make  the  run  and  reference  images  line  up.  The  run  image  will 
be  shifted  so  that  it  aligns  with  the  reference  image: 

I™„0'J')=I™„(i  +  Ai,j  +  Aj)  (4.1) 
where  Ai  and  Aj  are  the  amount  of  pixel  shifting  necessary.  Note  that  these  quantities  are 
generally  non-integer  values  in  which  case  a  weighted  average  of  the  four  closest  pixels 
will  be  used  to  calculate  the  required  value.  The  next  step  is  one  that  is  often  not  possible 
to  perform,  but,  for  the  experiment  to  be  presented  in  this  dissertation,  it  is.  A  tiny 
fluorescent  marker  was  applied  to  the  external  surface  of  the  test  section  in  a  location 
where  it  would  be  within  the  field  of  view  of  the  PSP  images.  Because  this  marker  is  on 
the  outside,  its  intensity  will  only  vary  if  the  intensity  of  illumination  varies.  Taking  a 
measure  of  the  average  intensity  over  a  region  of  the  marker  in  the  run  image  (minus  the 
average  dark  intensity  over  the  same  region)  and  ratioing  to  the  same  in  the  reference 
image,  the  run  image  can  be  scaled  to  match  the  intensity  of  illumination  that  was  present 
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during  the  acquisition  of  the  reference  image.  That  is,  define  the  "external  reference" 
scaling  constant: 

C.  =  |""'"|'"  (4.2) 

run.ext  drk 

When  multiplied  by  the  run  image,  this  constant  will  perform  the  necessary  scaling.  With 
constant  light  sources  such  as  the  tungsten-halogen  lamp  or  LED  lamp  used  in  the 
aforementioned  Dynamic  Test  Cell,  such  corrections  are  not  generally  necessary. 
However,  for  the  particular  application  to  be  presented,  a  strobe  lamp  is  being  used.  If 
the  shutter  on  the  camera  is  not  ultra  stable  in  its  opening  process,  closing  process  and 
total  duration  of  operation,  or  if  there  is  variability  in  the  strobe,  one  image  may  see  an 
additional  strobe  flash  or  two  relative  to  another  image  depending  on  the  fi-equency  of  the 
strobe  flashes.  So,  it  becomes  necessary  to  be  able  to  scale  out  the  intensity  of 
illumination  using  the  external  reference  technique.  It  is  being  assumed  here  that  the 
spatial  variations  due  to  a  single  flash  of  the  strobe  will  even  out  with  the  hundreds  of 
strobe  flashes  used  to  acquire  a  single  image.  The  next  step  is  to  filter  the  images 
spatially  to  remove  some  of  the  shot  noise  inherent  in  measuring  photonic  sources. 
Generally,  one  is  able  to  take  multiple  run  images  at  the  same  conditions  and  average 
them  to  reduce  the  shot  noise.  However,  for  the  particular  experiment  presented,  the  run- 
times were  limited  to  where  this  luxury  was  not  available.  Moreover,  the  integrated  light 
levels  of  many  strobe  flashes  was  still  not  sufficient  for  making  use  of  the  fiill  well 
capacity  of  the  CCD.  That  is,  the  measured  intensities  were  only  on  the  order  of 
hundreds  of  counts  (bits).  This  means  that  the  signal-to-noise  ratio  of  the  pixel  intensities 
is  relatively  low.  If  the  noisy  run  and  reference  images  were  ratioed,  the  noise  would 
only  be  amplified  on  those  pixels  that  just  happened  to  be  statistically  greater(lesser)  and 
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lesser(greater)  than  average  between  the  run  and  reference  images  respectively.  So,  at 
the  expense  of  some  spatial  resolution,  the  images  are  filtered  prior  to  being  ratioed.  The 
smoothing  kernel  used  was  of  the  form: 

cos'f7r^^7  +  j7/2nV2] 
K(io  Jo)=  ^;  ■  — —  where  i^,    =  -n,...,n;  n  =  1,2,3,...  (4.3) 

^  ^cos'(nVi'  +jV2nV2) 


i=-n  j=-n 


In  the  above  equation,  io  and  jo  are  the  coordinates  of  the  pixel  within  the  kernel  being 
calculated.  The  kernel  has  dimensions  (2n+l)  x  (2n+l),  thus  it  is  of  odd  order  (3x3,  5x5, 
. . .)  and  is  centered  on  the  pixel  ((n  + 1)/2 , (n  + 1)/2) .  It  may  be  noticed  that  the  comer 
pixels  of  this  kernel  will  turn  out  to  be  identically  zero.  This  kernel  is  then  2D- 
convoluted  with  an  image  to  achieve  a  smoothed  image.  This  form  of  kernel  is  superior 
to  a  3x3  or  5x5  "flat-top"  kernel  in  that  it  gives  more  weight  to  pixels  near  the  center  and 
less  weight  to  those  farther  away.  For  a  flat-top,  equal  weight  is  given  to  all  pixels  within 
the  area  of  support  of  the  kernel.  This  kernel  was  developed  based  on  ideas  found  in  ID- 
DSP  wherein  spectral  leakage  is  reduced  by  using  a  tapered  windowing  function  on  a 
dataset.  Depending  on  the  amount  of  noise  observed  (based  on  the  light  levels 
measured),  either  a  3x3  or  5x5  kernel  was  used  to  smooth  the  data.  An  estimate  for  the 
spatial  resolution  based  on  a  95%  sum  of  the  kernel  values  gives  a  value  roughly  equal  to 
"n".  As  an  example,  a  7x7  kernel  for  which  n=3,  roughly  95%  of  the  kernel's  weight  is 
contained  within  a  radius  of  3  pixels  about  the  center  of  the  kernel.  It  should  be  noted 
that  for  3x3  and  5x5  kernels,  100%  of  the  kernel's  weight  is  contained  within  the 
estimated  area.  The  next  step  in  the  flowchart  is  to  ratio  the  smoothed  images.  The  steps 
up  to  this  point  may  be  summarized  via  the  following  equations: 
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where  •  is  the  "array  multipHcation"  operator  ( C  =  A  •  B  ^     =  a^bij ),  and  (8i 

represents  the  convolution  operation.  Ratioing  these  two  equations  we  finally  end  up 
with  the  equation: 

■  _Y^_  K®[C..(I.,-1„)] 

Y„  -  K«[c„c„  .(I..  -i„r 

which  represents  the  static  calibration  function.  Analytically  inverting  this  function  gives 

Pru„=g(Y*,P„r)  (4.6) 
where  f  and  g  are  nearly  inverse  functions. 

One  fact  that  will  aid  in  the  analysis  of  the  images  for  this  particular  application  is 
that  we  have  the  capability  to  perform  a  static  calibration  of  the  PSP  sample  in  the  exact 
spot  and  with  the  same  lighting  as  during  a  run.  That  is,  the  exit  duct  may  be  unscrewed 
fi-om  the  test  section  and  replaced  with  a  cap.  The  test  section  is  now  an  enclosed  volume 
which  can  be  drawn  down  to  vacuum  or  can  be  pressurized.  Thus,  we  can  get  a  very 
accurate  measure  of  the  static  calibration.  Following  the  process  outlined  in  leading  up  to 
equation  (4.6),  a  sequence  of  images  taken  at  different  pressures  with  the  tunnel  capped 
may  be  reduced  to  a  sequence  of  data  points  of  dimensionless-intensity  ratio  versus 
dimensionless-pressure  ratio.  An  analytic  function  may  now  be  curvefit  to  this  data  to 
give  a  useful  expression  for  the  calibration  function.  Table  3  lists  mathematically 
equivalent  forms  of  three  common  calibration  functions.  The  first  fiinction  in  the  table  is 
the  standard  linear  Stem-Volmer  equation.  The  second  corresponds  to  a  linear 
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Table  3~Calibration  Functions,  Constraints,  and  Typical  Values 

Function  Constraints^  Typical  Values^''' 

.      .  (P-Q.)  Q,  <0.01  Q,  =-0.2 

°(P..-Q,) 


Y'-Y'(P-Q.XP.r-Q2) 
°(P.r-Q,XP-Qj 


Y*  =  y: 


Q,  <0.01 

Q,  =-0.1 

Q2<Q. 

Q,=-0.4 

Q3) 

Q,  <0.01 

Q,  =-0.05 

-Q3) 

Q2<Q. 

Q,  =  -0.3 

Q3<Q2 

Q3=-1.5 

Numerical  values  of  constraints  and  typical  values  are  given  in  units  of  atm. 
YJ  =  1.00  is  the  typical  value  for  all  three  functions. 

Stem-Volmer  response  with  a  pressure-independent  bias  in  intensity.  The  third  entry  is 
the  dual-mode  or  dual-sorption  calibration.  All  of  these  calibration  forms  have  the 
benefit  that  the  Qi's  will  be  independent  of  the  reference  pressure,  Pref.  The  calibration 
data  of  Y*  versus  P  for  the  PSPs  used  in  the  transverse  jet  experiments  had  a  non- 
negligible  amount  of  downward  curvature.  Therefore,  the  standard  Stem-Volmer 
calibration  was  insufficient  for  fitting  the  data.  It  was  found  that  the  third  function  from 
Table  3: 

Y'-Y-  (P-Q.XPrer-Q2XP-Qj 

"  "(P..-Q,XP-Q2XP„r-Q3)  ,  ^'-'^ 

most  accurately  fit  the  data.  The  parameter  Yj  is  only  used  during  the  curvefitting 
process.  It  will  be  a  value  very  near  1 .00  and  allows  for  uncertainty  in  the  reference 
image  during  the  curvefitting  process.  When  implementing  this  fiinction  to  calibrate 
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actual  PSP  data,  the  parameter  Yj  is  set  to  1 .  The  three  values  Qi,  Q2  and  Q3  must  all 
satisfy  certain  constraints  to  satisfy  physics,  otherwise  the  intensity  would  be  allowed  to 
be  infinite  at  some  small  positive  pressure.  The  three  constraints  listed  in  Table  3  are: 

•  Q,  <-0.01 

•  Q2<Q.  (4.8) 

•  Q3<Q2 

The  value  of -0.01  on  the  first  constraint  is  not  truly  a  physical  constraint  but  rather  an 
empirical  one.  Physically  speaking,  the  true  constraint  on  Qi  is  that  it  must  be  strictly 
less  than  zero.  However,  it  has  been  found  that  during  curvefitting,  this  parameter  will 
often  approach  the  upper  limit  and  perhaps  even  fit  to  the  upper  limit  value.  If  a  limit  of 
0.00  were  chosen,  the  calibration  would  be  indicating  that  when  the  pressure  goes  to  zero 
the  intensity  goes  to  infinity.  This  cannot  possibly  be,  so,  a  slightly  more  negative  value 
was  chosen  for  the  upper  limit  on  Qi.  One  additional  note  with  regard  to  these 
parameters  is  that  they  are  mathematically  the  poles  (Q2)  and  zeros  (Qi,  Q3)  of  the 
calibration  function.  The  inverse  of  equation  (4.7)  with  Yj  =  1  gives 

p  _  -B  +  Vb^  -4AC 
2A 

where 

A  =  (P„f-Qj  (4.9) 
B  =  -(P.,  -  Q, XQ,  +  Q3 ) -  Y- (P.,  -  Q.  Xp„,  -  Q3 ) 
C  =  Q,Q3(P„r  -Q2)+ Y-(P„,  -Q,Xp„,  -Q3)Q, 

A  sample  calibration  curve  is  shown  in  Figure  4.14.  The  data  spans  pressures  fi-om 
O.lSatm  up  to  1.67atm  (19kPa  -  ITOkPa).  The  curvefit  parameters  for  this  dataset  are 
(Qi,  Q2,  Q3,  Yj)  =  (-0.0100,  -0.  2254,  -1.  3211,  1.0066)  where  Q,,  Q2,  and  Q3  are  given 
in  units  of  atm.  Notice  that  Q,  has  gone  to  its  upper  constraint  and  that  Y*  is  very  close 
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to  1.  The  RMSE  of  the  curvefit  is  0.84%  of  P^r,  or,  0.85kPa.  If  the  constraints  are  lifted 
and  the  data  is  re-curvefit,  the  parameters  become  (0.0484,  -0.0759,  -1.0761,  1.0054)  and 
the  RMSE  is  0.77%  of  Pref,  or  0.78kPa.  So,  with  the  constraints  lifted,  the  curvefit  would 
indicate  that  the  intensity  is  infinity  when  the  pressure  is  P  =  0.0484atm ,  that  is,  when 
P  =  4.89kPa .  Notice  that  without  the  constraints,  the  curvefit  only  decreases  the  error  of 
fit  by  0.07kPa  which  is  a  very  negligible  amount  for  most  applications. 


For  a  function  f(x,y)  where  x  and  y  are  independent  variables,  the  general 
uncertainty  in  f  is  given  by 


For  the  independent  variables  x  and  y,  the  95%  confidence  limits  can  be  calculated  by  the 
formula 


where  B  represents  all  sources  of  systematic  or  bias  error  and  S  is  the  standard  deviation 
of  the  random  component  of  the  error.  Using  the  formulas  of  equations  (4.10)  and  (4.1 1), 
the  uncertainties  of  the  variables  for  the  pulsed  jet  experiment  can  be  tabulated.  Table  4 
lists  many  of  these  variables,  their  nominal  values  and  their  uncertainties. 

A  detailed  investigation  of  the  uncertainty  in  the  indicated  PSP  pressure  shows 
that  the  relative  uncertainty  in  the  measured  PSP  pressure  can  be  reduced  to  the  following 
equation: 


Uncertainty  Analysis 


(4.10) 


(4.11) 


(4.12) 


where  P  is  the  fiinction  p(y*  )  from  the  static  caHbration,  Y*  is  the  measured  intensity 
ratio,  P  is  the  A/D  gain  of  the  CCD  pixels  in  electrons  per  bit,  and  Iref  is  the  intensity  in 
the  reference  image.  A  full  derivation  of  this  equation  is  given  in  Appendix  C.  We  can 


Table  4~Nominal  values  and  95%  confidence  limits  of  key  variables. 


Variable 

Nominal  Value 

95%  Confidence  Limits 

P  stag 

1.555atm 

±0.021atm 

Pstat 

0.390atm 

+0.004atm 

Twall 

9.78°C 

±o.2rc 

M 

1.557 

±  0.006 

Pj 

2.000atm 

+0.014atm 

J 

1.134 

±0.057 

a 

13.5|jm 

±2.0^m 

E 

5.35kJ/mole 

±0.27kJ/mole 

T 

0.20s 

±0.06s 

Irun  "  Idfk 

260cts 

±  3.7cts 

Iref  -  Idrk 

160cts 

±  2.9cts 

Ppsp " 

0.390atm 

±0.016atm 

P 

0.390atm 

±  0.00  latm 

Assumes  no  smoothing  (see  Figure  4. 15); Upstream  and  downstream  pressure  taps. 
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further  take  into  account  the  effect  of  the  smoothing  kernel  (equation  4.3).  Smoothing 
has  the  effect  of  increasing  the  effective  number  of  counts  at  a  given  pixel  by  a  factor 
related  to  the  size  of  the  smoothing  kernel.  Specifically,  the  factor  is  given  by  the 
denominator  of  equation  (4.3).  That  is,  for  a  3  x  3  kernel,  the  effective  intensity  increases 
by  a  factor  of  2.032  and  by  a  factor  of  5.808  for  a  5  x  5  kernel.  The  result  of  this  in 
relation  to  the  relative  uncertainty  in  the  indicated  pressure  is  shown  in  Figure  4.15.  The 
effect  of  smoothing  in  relation  to  equation  (4.12)  is  also  given  in  Appendix  C. 


Figure  4.1--PSP  Dynamic  Test  Cell:  A)  Digital  valve  to  vacuum,  B)  Pressure 
transducer,  C)  Digital  valve  to  high  pressure,  D)  Thermistor,  E)  Window  to  PSP 
sample,  F)  Cooling  fins. 
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Figure  4.2~Schematic  view  of  the  PSP  Dynamic  Testing  Cell. 


Figure  4.3--Heat  Exchanger  Unit:  A)  Inlet,  B)  Outlet. 
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Figure  4.4--Oriel  Quartz-Tungsten-Halogen  (QTH)  Lamp. 


Figure  4.5--ISSI  LED  Lamp. 


PSP  Test  Cell 


Figure  4.6--Schematic  of  the  PSP  Dynamic  Testing  Cell  Optics. 


Figure  4.7--Fisher  Control  Valves:  A)  VA"  Valve,  B)  1"  Valve. 


Figure  4.8--Stagnation  chamber  of  the  supersonic  test  stand  facility. 


Figure  4.9-Mach  1.56  Test  Section:  A)  3  of  8  Static  pressure  taps,  B)  2  of  4 
Thermocouples,  C)  Exterior  reference  marker,  D)  Bottom  edge  of  circular  plug  insert. 


Figure  4.10~Detailed  front  view  of  plug  insert  showing  the  central  jet  orifice  and  the  4 
static  pressure  tap  locations. 


Figure  4.1 1 -Detailed  rear  view  of  plug  insert:  A)  2  of  4  Pressure  taps,  B) 
Thermocouple,  C)  Digital  valve  controlling  the  jet. 


Figure  4.12--A)  Xenon  Strobe,  B)  Photometries  CCD  Camera. 
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Figure  4.13--Image  Processing  Flowchart. 
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Figure  4.15--Relative  uncertainty  in  the  indicated  PSP  pressure. 


CHAPTER  5 
RESULTS 

Numerous  studies  have  been  done  on  the  effect  of  a  jet  of  air  injected  into  a 
supersonic  freestream/'^'"^"  All  of  these  tests  were  performed  using  a  steady  blowing  jet. 
As-yet,  there  have  been  no  experiments  involving  pulsed  jets  in  supersonic  crossflow. 
This  may  be  because  the  complexity  of  the  flowfield  under  steady  blowing  is  already 
difficult  enough  to  describe  and  quantify.  Or,  it  may  be  because  there  hasn't  been  a 
viable  measurement  technology.  However,  with  the  mounting  interest  in  synthetic  jets, 
which  on  small  time  scales  are  cyclic,  there  may  soon  be  great  interest  in  an  investigation 
of  a  pulsed  jet.  This  then  is  the  basic  idea  for  an  experiment  to  test  the  ability  of  the 
compensator  design  to  correct  wall  PSP  measurements  around  the  orifice  of  a  pulsed  jet. 

The  equipment  used  is  as  described  in  the  supersonic  wind  tunnel  section  of  the 
preceding  chapter.  The  tunnel  was  run  at  a  stagnation  pressure  of  158kPa  (23psi)  which 
was  the  lowest  value  possible  that  would  leave  the  tunnel  fully  started.  That  is,  if  the 
stagnation  pressure  dropped  much  below  that  value,  the  normal  shock  at  the  exit  of  the 
test  section  would  detach  from  the  exit  and  move  upstream  into  the  test  section  proper. 
The  reason  for  this  low  stagnation  pressure  was  to  allow  as  long  a  runtime  as  possible. 
Once  the  tunnel  had  been  started,  there  was  an  approximately  5  minute  equilibration 
period  during  which  the  wall  temperature  of  the  test  section  would  drop  over  time  and 
then  level  off  Depending  on  the  outside  temperature  (i.e.  the  temperature  of  the  exterior 
tanks)  the  static  wall  temperature  would  range  from  as  low  as  0°C  (in  the  wintertime)  up 
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to  more  normal  value  of  around  10°C  to  15°C.  It  should  be  noted  that  all  of  these 
experiments  were  performed  at  night.  This  was  because  ambient  light  from  outdoors 
would  come  through  the  exhaust  duct  and  into  the  test  section  creating  a  bias  intensity  in 
the  PSP  images. 

For  the  particular  tests  to  be  shown  here,  the  digital  valve  controlling  the  jet  was 
cycled  at  a  frequency  of  20Hz  with  a  40  percent  "open"  duty  cycle.  That  is,  the  jet  valve 
was  opened  for  20ms  and  then  closed  for  a  duration  of  30ms,  giving  a  total  period  of 
50ms.  The  strobe  could  be  set  to  any  desired  phase  relative  to  the  20Hz  jet  valve  control 
voltage.  A  single  PSP  image  had  a  duration  of  20  seconds  of  integration.  During  this  20 
seconds  the  strobe  would  have  flashed  a  total  of  20Hz  x  20s  =  400  flashes  at  exactly  one 
relative  phase.  Thus,  the  image  acquired  was  an  ensemble  average  of  400  flashes.  The 
strobe  would  then  be  reset  to  a  different  phase  and  another  image  would  be  taken. 

In  other  transverse  injectant  studies,  the  jet  pressiure  was  one  of  the  primary 
variables.  In  the  experiments  presented  herein,  the  jet  pressure  was  set  to  a  value  of 
around  2  atmospheres  absolute.  Based  on  this  value,  the  jet-to-crossflow  momentum 
ratio  (equation  1 .8)  was  J=l .  134.  After  a  full  sequence  of  images  had  been  taken,  the 
tunnel  was  quickly  shut  down.  As  soon  as  the  flow  in  the  tunnel  had  stopped,  a  reference 
image  at  ambient  pressure  was  acquired.  Because  this  reference  image  was  acquired 
soon  after  shutdown,  the  temperature  of  the  tunnel  wall  had  not  yet  had  time  to  rise  very 
far.  Typically,  this  temperature  rise  was  on  the  order  of  0. 1°C  to  0.2°C,  which  is  within 
the  limits  of  uncertainty  quoted  in  Chapter  4.  So,  a  reference  image  was  available  that 
was  taken  at  the  same  temperature  as  during  the  run  images.  Performing  the  calibration 
procedure  outlined  in  the  previous  chapter,  one  now  had  a  sequence  of  calibrated  PSP 
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images.  Because  of  the  very  low  light  levels,  a  5  x  5  smoothing  kernel  was  used  to  even 
out  the  pixel-to-pixel  statistical  variations.  Based  on  the  uncertainty  analysis  presented  in 
Chapter  4  and  in  Appendix  C,  the  uncertainty  of  the  PSP  indicated  pressures  will  be 

reduced  by  a  factor  of  V5.808  from  the  value  presented  in  Table  4.  From  that  point,  the 
experiment  now  proceeds  into  the  compensation  phase. 

The  time  constant  of  the  PSP  was  calculated  from  a  compilation  of  various 

experiments.  Recalling  the  definition  of  the  time  constant  from  Chapter  2,  x  s  a      ,  and 
assuming  an  Arrhenius  relation  for  the  diffiisivity,  equation  (1.12),  the  necessary 
measurements  are  evident.  The  thickness  was  measured  at  20  points  over  the  surface  of 
the  plug  insert  using  an  eddy-current  gauge  (Positector  6000)  and  then  averaged.  At  the 
same  time  the  plug  insert  was  painted,  two  small  samples  were  also  painted  for  use  in  the 
Dynamic  Test  Cell  presented  in  Chapter  4.  The  PSP  on  these  samples  were  roughly  0.4x 
and  1 .2x  the  thickness  of  the  coating  on  the  plug  insert.  The  thickness  of  these  coatings 
was  also  measured  using  the  eddy-current  gauge.  Using  the  Dynamic  Test  Cell,  a  series 
of  tests  were  performed  on  the  thicker  of  the  two  small  samples  to  measure  the  time 
constant  of  the  PSP  as  a  fimction  of  temperature.  The  data  was  then  curvefit  using  the 
formula 

T,  =Toexp(E/RT),  (5.1) 
where  E  is  the  activation  energy  of  diffiision.  Because  the  thickness  of  this  coating  was 
not  equal  to  the  coating  thickness  on  the  plug  insert,  a  fiirther  step  was  needed  to 
calculate  the  value  of  the  time  constant  to  be  used  in  the  compensation  routine.  The 
formula 


X  = 


a 
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X,  (5.2) 


corrects  the  time  constant  for  differences  in  coating  thickness  between  the  small  sample, 
a, ,  and  the  plug  insert,  a.  Using  the  measured  wall  temperature  in  combination  with 
equations  (5.1)  and  (5.2),  a  value  of  x  =  0.20s  was  calculated  for  the  PSP  on  the  plug 
insert  and  used  for  the  purposes  of  dynamic  compensation. 

Shown  in  Figure  5.1  are  a  sequence  of  10  calibrated,  but  uncompensated,  PSP 
images  taken  at  equal  time  steps  throughout  the  cycle.  In  terms  of  a  temporal  measure, 
the  images  were  separated  in  time  by  5ms  per  image.  The  image  in  the  upper  left-hand 
comer  of  the  figure  corresponds  to  the  opening  of  the  jet  valve.  Progressing  from  the  top 
down  and  then  to  the  top  of  the  second  column  and  down  spells  out  a  fiill  cycle.  The 
phase  delays  relative  to  the  opening  of  the  jet  valve  are  shown  on  each  image  as  a 
percentage  of  a  full  cycle.  The  gray-scale  of  the  images  has  been  set  so  that  white 
corresponds  to  a  pressure  of  0.7atm  and  black  is  0.1  atm.  It  can  be  seen  that  the  area 
upstream  of  the  jet  (left  of  the  jet)  is  generally  a  region  of  high  pressure  and  that  the  area 
downstream  of  the  jet  (to  the  right)  is  a  region  of  low  pressure.  From  the  sequence  of  10 
images,  the  upstream  pressure  indicates  a  peak  pressure  of  between  0.55  and  0.60atm  for 
phase  delays  of  10%  to  50%.  The  minimum  pressure  downstream  of  the  jet  is  between 
0.25  and  0.28atm  and  occurs  around  a  phase  delay  of  50%.  The  mean  pressure  far  from 
the  jet  is  on  the  order  of  0.39atm.  So,  based  on  the  uncompensated  PSP  images,  the  peak 
wall  pressure  is  0.19atm  above  the  freestream  and  the  minimum  wall  pressure  is  0.13atm 
below  the  freestream.  Note  that  the  uncertainty  in  these  values  is  approximately 
0.007atm.  So,  the  uncertainty  is  on  the  order  of  5%  of  the  absolute  fluctuation. 
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Progressing  now  to  the  compensation  data,  the  sequence  of  10  compensated 
images  is  shown  in  Figure  5.2.  These  images  were  generated  using  the  Matlab  and 
Simuhnk  codes  presented  in  Appendices  A  and  B.  Running  on  a  Pentium  III  800MHz 
machine  with  256MB  of  RAM,  compensating  a  sequence  of  10  images  (177pixels  x 
253pixels)  took  approximately  2hrsl5min  to  compile.  In  these  compensated  images  we 
see  the  shock  strength  much  clearer.  Also  note  that  the  pressure  very  quickly  dies  down 
after  the  jet  is  shut  off  as  compared  to  the  uncompensated  where  it  appears  the 
downstream  pressure  never  equilibrates.  The  pressure  midway  between  the  front  edge  of 
the  jet  orifice  and  the  upstream  pressure  tap  peaks  at  around  0.70atm  at  a  phase  delay  of 
0%.  The  pressure  downstream  of  the  jet  reaches  a  minimum  around  0.1  atm  almost 
immediately  downstream  of  the  jet  at  a  phase  delay  of  20%.  This  gives  a  peak  wall 
pressure  0.31  atm  above  the  freestream  and  a  minimum  of  0.29atm  below  freestream. 
Another  interesting  feature  is  that  at  the  initiation  of  the  jet  (0%  phase  delay)  the  whole  of 
the  image,  both  upstream  and  downstream,  is  seen  to  raise  in  pressure.  This  seems  to 
indicate  that  in  the  initial  instant  of  the  jet,  it  is  filling  up  the  boundary  layer  in  all 
directions,  pumping  it  up  to  a  higher  pressure.  Then,  the  gross  flow  is  able  to  overcome 
this  initial  transient  and  set  up  a  more  clearly  defined  shock  structure.  This  shock 
structure  persists  nearly  unchanged  until  the  jet  closes.  When  the  jet  is  then  turned  off, 
the  flow  nearly  immediately  adjusts,  wiping  out  all  evidence  of  pressure  gradients. 

Shown  in  Figure  5.3  is  a  comparison  of  the  pressures  along  the  streamwise 
centerline  through  the  jet  for  both  the  uncompensated  and  the  compensated  images.  The 
centerline  pressures  from  a  set  of  10  images  were  stacked  into  a  2D  array.  This  2D  array 
was  then  duplicated  and  stacked  twice  yielding  a  three-cycle  loop,  thus  giving  a  clearer 
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Figure  5.1--A  sequence  of  10  uncompensated  PSP  images  taken  at  5  ms  intervals  during 
a  full  cycle  of  50ms.  Grayscale  is  such  that  white  =  0.7atm  and  black  =  0.  latm. 
Dimensions  are  from  ^.5Dj  to  +5.5Dj  in  X  and  from  -3.5Dj  to  +3.5Dj  in  Y. 


Figure  5.2--A  sequence  of  10  compensated  PSP  images  taken  at  5  ms  intervals  during  a 
full  cycle  of  50ms.  Grayscale  is  such  that  white  =  0.7atm  and  black  =  0.  latm. 
Dimensions  are  from  -4.5Dj  to  +5.5Dj  in  X  and  from  -3.5Dj  to  +3.5Dj  in  Y. 
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Figure  5.3--A  three-cycle  loop  of  PSP  data  along  the  streamwise  centerline  through  the 
jet.  (a)  Uncompensated;  (b)  Compensated.  Grayscale  is  such  that  white  =  0.7atm  and 
black  =  0.  latm.  The  X-coordinate  is  oriented  along  the  horizontal  direction  (positive  to 
the  right)  and  time  is  oriented  vertically  (positive  up). 
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Figure  5.4--Comparison  of  PSP  data  to  pressure  tap  data,  (a)  Upstream  of  the  jet;  (b) 
Downstream  of  the  jet;  (c)  Spanwise  from  the  jet. 
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Figure  5.4--continued. 


Figure  5. 5-- Variation  in  the  paint  thickness  near  abrupt  edges. 
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visual  presentation  of  the  data.  The  enhanced  features  due  to  the  compensator  are  quite 
evident. 

Shown  in  Figure  5.4  are  comparisons  of  pressures  measured  by  the  pressure  taps 
to  the  measurements  of  the  PSP,  both  compensated  and  uncompensated.  Once  again,  a 
single  cycle  has  been  looped  3  times  to  give  a  clearer  presentation.  The  pressure  tap  data 
is  an  ensemble  average  of  10  cycles.  The  PSP  data  is  an  average  of  32  values  taken 
around  the  pressure  tap.  The  32  values  are  at  equal  angular  spacing  at  a  radius  of  14 
pixels  (1.6mm;  0.63Dj)  from  the  center  of  the  tap.  This  large  radius  was  necessary  to  get 
away  from  the  corrupted  data  local  to  the  tap  as  will  be  explained  shortly.  In  the  figure,  a 
large  amount  of  overshoot  is  seen  at  the  taps  just  upstream  of  the  jet  and  spanwise  from 
the  jet.  This  is  likely  due  to  an  overestimation  of  the  time  constant  of  the  paint.  That  is, 
increasing  the  time  constant  used  in  the  compensator  relative  to  the  "true"  value  of  the 
time  constant  tends  to  cause  instability  in  the  compensator  and  yields  data  with  both 
overshoot  and  ringing.  Underestimating  the  time  constant  does  not  have  this  effect  and  is 
thus  a  preferable  situation,  however,  it  too  will  yield  an  incorrect  amount  of 
compensation.  Another  possible  explanation  for  the  overshoot  is  that  we  only  have  10 
samples  per  cycle.  During  compensation,  the  Simulink  routine  of  Appendix  B  runs  on  a 
finer  temporal  grid  spacing  by  linearly  interpolating  between  the  data  it  is  given.  If  there 
is  large  uncertainty  in  the  data  given  to  the  Simulink  routine,  which  is  true  in  this  case, 
there  can  be  large  discontinuities  in  the  time  derivative  of  the  data.  The  compensator  will 
then  amplify  these  resulting  in  data  that  appears  to  have  overshoot  or  ringing. 

The  data  immediately  surrounding  the  jet  is  corrupted  for  one  primary  reason  - 
paint  thickness  variations.  Figure  5.5  shows  an  image  of  the  PSP  surface  taken  at  a 


glancing  angle.  It  can  be  seen  that  the  paint  thickness  around  the  taps  and  the  jet  orifice 
varies  to  a  noticeable  degree.  Starting  from  the  edge  of  a  tap,  the  paint  increases  in 
thickness  as  the  distance  from  the  tap  increases.  Then,  at  a  point  approximately  7  pixels 
(0.8mm;  0.27Dj)  fi-om  the  center  of  the  tap,  the  paint  reaches  a  maximum  thickness.  The 
thickness  then  decreases  slightly  down  to  the  nominal  thickness.  These  thickness 
variations  are  due  to  surface  tension  effects  at  the  time  of  paint  application.  Because  the 
time  constant  of  the  paint  is  a  strong  function  of  the  thickness,  the  compensated  data  in 
the  extreme  near-field  region  of  the  taps  cannot  be  trusted. 

The  last  comparison  to  be  done  is  that  of  the  effective  back-pressure  ratio  shown 
in  Figure  5.6.  This  data  was  calculated  in  a  similar  manner  to  the  data  surrounding  the 
pressure  taps.  At  a  radius  of  19pixels  (2.2mm;  0.74Dj)  fi"om  the  jet  center,  32  values 
around  the  circumference  of  the  jet  orifice  were  averaged.  This  value  was  then  divided 
by  the  pressure  behind  a  normal  shock  in  the  Mach  1.56  fi^eestream.  The  calculated  ratio 
shown  in  the  figure  is  nominally  on  the  order  of  0.39  except  for  a  10ms  window 
immediately  upon  the  opening  of  the  jet  at  which  time  the  ratio  peaks  at  around  0.50.  It 
should  be  noted  that  the  value  of  J=  1.1 34  quoted  earlier  is  only  valid  during  the  period 
when  the  jet  is  open.  When  the  jet  is  shut  off,  the  momentum  ratio  should  go  to  zero  and 
the  effective  back-pressure  ratio  should  approach  a  value  of  0.379.  This  back-pressure 
ratio  would  be  a  characteristic  value  for  the  M=1.56  flow.  Evaluating  the  uncertainty  in 
the  measured  back-pressure  ratio  reveals  an  uncertainty  of  0.017  relative  to  the  value  of 
0.39.  Thus,  the  characteristic  value  of  0.379  is  within  the  level  of  uncertainty  present  in 
the  experiment.  A  detailed  presentation  of  the  effective  back-pressure  ratio  as  a  function 
of  the  momentum  ratio  was  performed  by  Everett  et  al.*  Fitting  a  straight  line  to  their  data 


and  extrapolating  to  a  value  of  J=  1.1 34  yields  a  value  for  the  effective  back-pressure  o 
0.339.  This  value  does  not  compare  particularly  well  to  the  values  measured  here, 


however,  it  does  compare  well  in  relative  to  values  from  other  studies 
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Figure  5.6~Effective  back-pressure  ratio  at  the  jet. 
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CHAPTER  6 
CONCLUSIONS  AND  RECOMMENDATIONS 

Conclusions 

A  thorough  investigation  of  the  response  of  PSP  to  unsteady  pressures  was 
presented.  A  physically  realistic  model  for  the  dynamic  behavior  of  the  coatings  was 
developed  and  then  used  to  design  a  compensation  scheme.  This  compensation  scheme 
was  based  on  the  Model  Reference  Controller  and  was  in  effect  designed  to  be  state 
observer.  Compensation  of  data  acquired  under  controlled  laboratory  conditions  yielded 
excellent  results.  The  ambitious  task  of  applying  the  compensation  scheme  to  a  sequence 
of  full-field  images  was  undertaken.  Owing  to  the  large  experimental  uncertainties, 
particularly  in  the  overall  intensity  from  the  paints  and  from  the  spatially  varying 
thickness  of  the  PSP  coating,  the  compensated  images  could  only  be  used  in  a  semi- 
quantitative sense.  However,  the  designed  compensator  has  the  capability  to  be  a  great 
advancement  in  the  use  of  PSPs. 

Recommendations 
From  the  uncertainty  analysis  performed  for  the  full-field  dynamic  PSP 
experiment  presented,  the  following  recommendations  are  suggested  to  increase  the 
accuracy  of  the  PSP  measurements. 

1 .  Re-design  the  existing  test  section  to  reduce  it  in  size.  This  will  have  the  effect  of 
lowering  the  mass  flowrate  through  the  tunnel.  If  the  correct  size  is  chosen,  it  will 
allow  the  existing  compressor  system  to  keep  up  with  the  demands  of  the  tunnel 
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and  result  in  the  ability  to  run  the  tunnel  continuously.  With  unlimited  runtimes, 
many  more  cycles  of  the  pulsed  jet  can  be  summed  and  averaged,  thus  giving  the 
PSP  images  a  much  improved  signal-to-noise  ratio. 

2.  Develop  a  method  to  measure  the  coating  thickness  on  a  pixel-by-pixel  basis. 
This  will  result  in  a  better  estimate  for  the  time  constant  of  the  paint.  From  the 
current  experiments,  a  single  value  of  the  time  constant  was  used  for  the  entire 
image.  Allowing  the  time  constant  to  vary  spatially  should  result  in  much  better 
compensations. 

3.  Install  a  pressure  tap  inside  the  jet  between  the  valve  and  the  jet  orifice.  Then, 
install  another  Kulite  transducer  at  this  location  to  get  a  better  measurement  for 
the  jet  pressure.  This  should  then  allow  a  calculation  to  be  made  as  to  whether  the 
jet  is  sonic  at  the  exit  at  every  point  in  the  phase  of  the  jet. 

4.  Look  for  a  brighter  strobe  or  perhaps  one  with  a  slightly  longer  flash  duration. 
Currently,  the  flash  duration  is  about  an  order  of  magnitude  shorter  than  it  needs 
to  be  to  freeze  the  flow  for  the  range  of  fi-equencies  capable  by  the  jet  valve. 

5.  Acquire  more  than  just  10  images  per  cycle.  Finer  temporal  resolution  will  result 
in  less  noise  in  the  compensation  process.  A  minimum  of  20  images  per  cycle  is 
recommended. 

6.  A  final  recommendation  is  to  rotate  the  plug  by  a  small  angle  of  perhaps  10  or  15 
degrees.  Assuming  the  flowfield  is  symmetric  about  the  streamwise  centerline, 
the  paint  data  from  one  side  of  the  centerline  can  be  used  to  approximate  the 
pressure  indicated  by  the  tap  on  the  opposite  side  of  the  centerline.  Thus 
eliminating  the  need  for  averaging  the  paint  data  around  a  tap. 


APPENDIX  A 
MATLAB  CODE:  PSPMRCPARAMS 

%  PSP  Linear  and  Nonlinear  System  Models 
%  and  Model  Reference  Controller  (MRC) 
%  for  Importing  into  Simulink 

format  short  e 


%  basic  parameters 

%   

ord=5;  %  SS  Model  Order 

nx=20;  %  #  X  Grids  minus  1   (must  be  even) 

dx=l/nx;   %  X  Grid  Spacing 

tm=0.01;  %  Model  System  Time  Constant 

Q. 
"S 

%  calc  lin  and  nlin  systems  SS  form 
%   

a=zeros (ord) ; 
b=zeros  (ord,  2 )  ; 
c=zeros (nx+1, ord) ; 
d=zeros (nx+1 , 2 ) ; 
alin=zeros (ord) ; 
blin=zeros (ord, 2 ) ; 
clin=zeros ( 1, ord) ; 
dlin=zeros (1,2) ; 
for  i=l:ord, 

lam(i)=(2*i-l) *pi/2; 

a (i, i) =-lam(i) ^2; 

alin (i, i) =a  (i, i)  ; 

b(i,2)=-2*sin(lam(i) ) /lam(i) ; 

blin(i,2)=b(i,2)  ; 

clin(l,i)=-b(i,2) /2; 

for  j=l:nx+l, 

x=(nx+0.999999-j) *dx; 

%  use  0.999  instead  of  1  so  that  xmax<l 

c ( j , i) =cos (lam(i) *x) ; 

d(j,l)=l; 
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end 
end 

dlin=[l,0] ; 


%  convert  from  SS  form  to  TF  form 

o 
o 

[numl, denl ]=ss2tf(a,b,c,d,l) ; 
[num2,den2]=ss2tf (a,b,c,d,2) ; 

[numlinl, denlinl] =ss2tf (alin, blin, clin, dlin, 1 ) 

[numlin2, denlin2 ] =ss2tf (alin, blin, clin,  dlin, 2 ) 

%  num2(:,l)==0  so  remove  them 

len=size (num2 , 2 ) ; 

num2=num2 ( : , 2 : len) ; 

%  numlin2 (1) ==0  so  remove  it 

len=length (numlin2) ; 

numlin2=numlin2 (2 : len)  ; 


%  convert  from  MIxO  w/  U  =  [P,Pdot] 

%  to  SIxO  w/  U  =   [P],   i.e.  get  rid 

%  of  Pdot(t)   as  an  input 

%  

for  i=l:nx+l, 

num ( i , : ) =den2+conv ( [ 1  0 ] , num2 ( i , : ) ) ; 

end 

den=den2 ; 

numlin=denlin2+conv ( [1  0] , numlin2) ; 
denlin=d  :i2; 

"5 

%  define  linear  part  of  nlin  plant 


[ z ,  p ,  k ] =t  f 2  zp ( num, den ) ; 

%  ;  

%  define  lin  plant 

%  

[zp, pp,  kp] =tf 2zp (numl in, denl in) ; 


%  define  augmented  system  using  next 

%  highest  pole  in  sequence  from  denlin(s) 


%  note:  0(1)   system  =  ka/(s+ka) 

Q, 
O 

ka=(  (2*  (ord+l)  -1)  *pi/2)  ^^2; 
numa= [ ka] ; 
dena= [1  ka] ; 

[ za, pa,  ka] =tf 2zp (numa, dena) ; 

o 
o 

%  augment  lin  plant  w/  fast  0(1)  system 

o 
o 

numlin=conv (numa, numlin) ; 
denlin=conv (dena, denlin) ; 


%  make  numlin (s)   a  monic  polynomial 
%  

klin=numlin (1) ; 
numlin=numlin/klin; 
dp=length (denlin) -1 ; 

%  (klin,  numlin,  denlin)  define  augmented 
%  plant  for  MRC  calculation  purposes 

0, 

G   

■  define  model  system 

a 

km=l/tm;  '. 
numm= [ 1 ] ; 

denm=  [  1  km]  ;  -  , . 

Q, 

"O   — — 

%  dO  is  of  degree  dp  and  is 

%  arbitrarily  chosen  as  fast  or  faster 

%  than  the  model  system 

%  

dO=l; 

for  i=l:dp 

dO=conv([l  10* (i+9) ] ,dO) ; 

end 

o 
o 

%  calc  matrix  aa 


%  

%  column  1  of  aa 

acl= [denlin, zeros (1, dp-1) ] ; 

%  column  dp+1  of  aa 

ac2= [0, numlin, zeros (1, dp-1) ] ; 

for  1=1: dp, 

aa (i, : ) =acl ; 

acl= [0, acl] ; 

acl=acl  (1:  (2*dp) ) ; 

aa (i+dp, : ) =ac2; 

ac2=[0,ac2] ; 

ac2=ac2 (1: (2*dp) ) ; 

end 

aa=aa '  ; 

%  

%  calc  vector  bb 

%   

bb=conv (numm, denlin) -conv (numlin, denm) 
bb=conv (dO, bb) ; 

%  note  that  bb(l)==0  so  remove  it 
len=length (bb) ; 
bb=bb(2:len) ; 
bb=bb ' ; 

O       —  —  —  —  —  —   —  —  —  —  —  —  —  

%  calc  soln  vector  xx 

%  

xx=aa\bb; 

relerr= (aa*xx-bb) . /bb 

%  

%  calc  feedback  vectors  cl  and  c2 
%  

cO=km/klin; 
cl=xx ( 1 : dp) ; 
c2=xx(dp+l:2*dp)  ; 
[zl,pl,kl]=tf2zp(cl' ,dO)  ; 
[z2,p2,k2]=tf2zp(c2',d0) ; 


%  Simpson's  rule  vector  coeffs 
%  integrated  intensity  = 


%  (simp)  dot  (int  of  x) 
%  


for  i=l:nx+l, 

if  mod(i,2)==0 

simp (i) =4 ; 
else 

simp (i) =2; 

end 

end 

simp (1) =1; 
simp (nx+1) =l; 
simp=simp/sum  (simp)  ; 

o 
o 

%  calc  butterworth  filter  params 

a 


[zb,pb,)cb]=butter(5,l/tm,  's' )  ; 


APPENDIX  B 
SIMULINK  CODE:  PSPMRC 


Main  Routine:  pspmrc.mdl 


o- 


CD- 


f(u) 

1  - exp ( ■ I "  2 


PToducl4 


te-za(s) 
pa(8) 


PIX«)         ^  "J 


l«'za(i) 
*  P»(t) 
■ugL.iyilA'nl 


k2'z2|s) 
p2(l) 


-f5- 


PiDductl 


Subroutine:  nlin  system.mdl 


CD-H 

p 


02  iy«em 

(1-d  diffuaon) 


p2int 


coefft 


Dot  Product 


pv 


Subroutine:  p2int.mdl 


CD— +Q — 
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•    Subroutine:  int2p.mdl 


•  Notes: 

1 .  The  following  Simulink  blocks  are  continuous  time  systems  in  pole-zero  form: 
aug  system,  HI,  H2,  butterworth,  02_system. 

2.  The  following  parameters  are  loaded  from  the  Matlab  workspace  into  Simulink:  ka, 
za,  pa,  kp,  zp,  pp,  kl,  zl,  pi,  k2,  z2,  p2,  km,  k,  z,  p,  q,  simp. 

3.  Follow  these  steps  from  within  the  Main  Routine,  pspmrc.mdl,  in  order  to  load  the 
parameters  from  Note  2:  Simulation  >  Parameters  >  Real-Time  Workshop  >  Inline 
Parameters  >  Tunable  Parameters.  Enter  each  of  the  variables  from  Note  2  as  a 
tunable  parameter. 


APPENDIX  C 
PSP  UNCERTAINTY  ANALYSIS 

What  follows  is  a  derivation  of  equation  (4. 12).  It  should  be  noted  that  the  result 
derived  here  is  only  valid  for  the  case  of  low  light  output  from  the  PSP.  That  is,  when  the 
light  levels  are  low  the  uncertainty  due  to  photon  counting  statistics  will  be  the  dominant 
source  of  error  in  the  measurement.  All  other  sources  of  uncertainty  will  be  ignored  as 
they  are  assumed  to  be  negligible  in  comparison  to  the  photonic  noise.  To  begin, 
consider  the  uncertainty  in  the  indicated  PSP  pressure 


Up  = 


fdp  ^ 

2 

(d?  ] 

+ 

— u, 

Ui  'J 

1 51.  'J 

(C-1) 


where  I  is  the  intensity  of  the  run  image  and  Ir  is  the  intensity  of  the  reference  image.  If 
we  define  Y  =  I^/I,  the  chain  rule  can  be  employed  to  change  the  partial  derivatives  in 
equation  (C-1)  to  arrive  at 


Up  = 


ap  I 

5YI 


7^. 


+ 


dP  1 


-u, 


aYi 


(C-2) 


Simplifying  (C-2)  yields 


u  =Y 
p 


dY] 


hi 

2 

+ 

uJ 

J 

(C-3) 


The  full-well  capacity  of  the  pixels  of  the  CCD  camera  used  in  the  experiments  presented 
herein  is  320,000  electrons.  The  CCD  further  has  a  14-bit  digitizer  on  the  A/D  process. 
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Dividing  the  total  electrons  by  the  number  of  discrete  levels  of  the  AID  gives  a  pixel  gain 
of  the  CCD  of 

p  =  19.53  e'/bit.  (C-4) 
The  intensities  measured  by  the  CCD  camera,  i.e.  the  counts  or  bits,  is  actually  then  a 
number  of  electrons  counted 

e"  =  pi .  (C-5) 
Photonic  sources  may  be  described  by  the  Poisson  distribution.'*^  For  Poisson  statistics, 
the  variance  is  equal  to  the  square  root  of  the  measurement.  Therefore,  the  uncertainty  in 
the  number  of  electrons  is  given  by 

u,=V7  =  VpI  (C-6) 
which  gives  the  uncertainty  in  the  intensity  as 


1^=  1 


(C-7) 


Beginning  from  equation  (C-3)  we  can  then  say  the  following 


Un    =  Y  , 

"  BY 


ay 


dY 


I 


+ 

1 1. 

J 

J 

(  rrV  ( 
+ 


1 


V 


fx 

PIr 


+ 


) 


=  Y 


1  +  Y 

dYi  pi, 


(C-8) 
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Or,  we  have  the  relative  uncertainty  as 


Up  _  Y  ap  [l+Y 


(C-9) 


P     ?  dY\  pi. 


Equation  (C-9)  is  identically  the  equation  (4.12)  which  ends  the  proof. 

One  additional  point  that  can  be  added  is  the  effect  of  smoothing  the  images.  As 
mentioned  previously  in  Chapter  4,  the  smoothing  kernel  has  the  effect  of  boosting  the 
effective  number  of  counts  at  a  give  pixel,  thus  enhancing  the  statistics.  We  will  define  k 
as  the  denominator  of  equation  (4.3), 


The  variable  Y  will  be  unaffected  as  it  will  multiply  into  both  the  numerator  and  the 
denominator.  However,  it  will  scale  the  denominator  of  the  square  root  of  equation  (C-9) 
yielding 


Thus,  the  smoothing  kernel  will  improve  the  uncertainty  in  the  measured  pressure  by  a 


(C-10) 


(C-U) 


factor  of  Vic . 
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